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JAN BLATON 
1907—1948 


Obituary notice by A. RUBINOWICZ, Warsaw. 


7 The tragic death of Professor Jan Bla- 
ton on 17th May 1948 as the result of an 
accident in the Tatra Mountains was a severe 
blow to all of us. We have lost in him, 
a physicist, devoted heart and soul to scien- 
tific research, a teacher who carried away 
his students and who kindled in them the 
flame of desire for independent scientific 
thought, a man with a character like rock 
erystal. The knowledge that this eruel blow 
need never have fallen and that a little: 
more presence of mind might have sufficed 
to prevent the terrible catastrophe has 
alarmed us and has, if possible, plunged 
us all into even greater sorrow. 

Jan Blaton was born on 17th May 1907 in Sporysz in the beautiful‘ 
hills near the locality of Żywiec. He was.the fifth and youngest child 
-of a worker’s family. Presumably it never entered his head that the 
mountains in which he had spent his childhood, those mountains 
which he loved so dearly and which in his young days he so passionately 

plored, were plotting his untimely death. 

After his having graduated, with honours, at the gymnasium 
i Bielsko, we find Blaton enrolled in 1925 at the Faculty of Civil 
E::gineering at the Lwów Institute of Technology. He was, however, 
more interested in theoretical problems than in technical questions 

and therefore he passed to the General Faculty of this Institute in order 
to study physies. It was here, in 1928, that I made his aequaintance. 
. He was an uncommonly capable student, and his progress in the study 
of theoretical physies was astonishing. Even so, I was. quite taken 
aback when this young enthusiast for theoretical physics came before 
the summer holidays after barely one year of studies to ask me for 
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a subject for this magister examination thesis. I was, however, more 
astounded still when, after the holidays, he showed me the thesis 
almost completed. This work, on the intensities of the quadrupole 
lines in the Balmer series, was the first of a series of his publications 
concerning multipole radiation, among which we find his most valuable 
scientific achievements. 

In 1929 he became my assistant. He held this position only for 
a few weeks, however, since he was at that time under court sentence 
for two weeks detention —the sentence remaining suspended — for 
the distribution of communist leaflets. Devoted, to the point of self- 
sacrifice, to the cause of the progress of humanity, Blaton had entered 
directly after he had passed his matriculation into contact with 
a group of students belonging to the students’ union „Życie“ (Life) 
and in spite of intense scientifie work he continued to be politically 
active. Considerable efforts were needed to make it possible for him 
to remain at the Institute, if only in the character of a student. 

In March 1931 Blaton took his magister degree and less than two 
months later he submitted a second paper for publication. Its purpose 
was to investigate the problem of whether, in the scattering of light 
quanta by atoms, two photons can be turned into one photon. 

Blaton’s third paper was presented to the Polish Academy of 
Sciences and Letters in July of the same year. In it the author en- 
deavoured to prove that the electric field of a wave of light, despite 
its periodicity, causes a Stark effect. 

In 1932, Blaton was awarded a National Culture Fund Scholarship 
and, on the basis of his thesis on the dispersion of light in the vinieity 
of quadrupole lines, he took his doctor's degree at the Lwów Institute 
of Technology. Together, in this same year, he and I wrote a report 
on quadrupole radiation for the ,Ergebnisse der exakten Natur- 
wissenschaften“. __ 

Continuing to benefit from the National Culture Fund Scholarship, 
Blaton proceeded to Munich where, under the direction of Professor 
Sommerfeld, he continued his studies until the advent of Hitlerism 
forced him to leave Germany. He then went to Zurich. 


From 1933 to 1935, after his return to his native land, Blaton — 


occupied the post of an assistant in the Department of Theoretical 
Physics at the University of Wilno. One of his best scientifie achieve- 
ments, that is his participation in the discovery of magnetic dipole 


lines, dates from this period. In 1933 H. Niewodniczariski discovered _ 


that in mixtures of lead vapour with helium and argon there appear 
in the spectrum of high frequency electric discharges forbidden lines 
of non-ionized lead. In a letter dated November 1933, addressed to 
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the Editors of the Physical Review, Blaton together with Niewodni- 
czański came to the conclusion that one of these lines must be a magnetic 
dipole line not containing any admixture of electric quadrupole radi- 
ation. Niewödniezanski, by investigating its Zeeman effect, confirmed this 
result later and thus finally discovered the magnetic dipole radiation. 

The year 1934 saw Blaton submitting his thesis for admission 
as lecturer of Theoretical Physies at the Wilno University. In it the 
author shows how to calculate the intensities of magnetic dipole lines in 
case of any coupling. 

In 1935 Blaton published a paper on quaternions, semivectors 
and spinors, investigating the relations between these quantities. 

Blaton was appointed, in 1936, Head of the National Meteorol- 
ogical Institute. He occupied this post until the outbreak of the war 
A proof to the versatility of his intellect is the rapidity with which 
he mastered meteorology, a sphere of science at that time quite unknown 
to him. Two papers of his, quoted in meteorological textbooks, provide 
evidence of this. In one, the author investigates the dependence of 
the wave lengths of the gravitational waves on depth, employing 
methods used in wave mechanics. In the other, he gives for any .plane 
fluid motion a simple kinematic relation between the rate of change 
of the velocity direction of the fluid motion and the radii of curvatures 
of the path of the fluid particles and of the stream line. Applying 
this formula to the dynamics of the atmosphere, he demonstrates 
how it is possible to calculate changes in the direction of winds with 
a horizontal movement in the atmosphere with a knowledge only 
of the coefficient of friction and the data derived from synoptic maps. 

Since the post he then held involved a considerable degree of 
responsibility, it was necessary for Blaton to concentrate his whole 
mind and attention on meteorology; yet he found the time, in 1937, to 
write the last of his multipole papers. In it he explains how to split 
the radiation created by any periodic electric currents into electric 
and magnetic multipole radiations. Unfortunately, this paper appeared 
only in the Acta Physica Polonica and was overlooked by certain 
foreign authors. 

During the German occupation of Poland, Blaton ou at the 
Jegiel District Forestry House and, in order to lecture on Theoretical 
Mechanies and Physies at the clandestina University, journeyed to 
Warsaw at regular intervals until July 1944. His enthusiasm in Over- 
coming difficulties and the spontaneous admiration he aroused in 
his students are proved by the fact that amidst the difficult war con- 
ditions he managed to educate several pupils who have have already. 
become independent scientific workers. 
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During the war years Blaton was also occupied with writing 
a textbook of mechanies. Judging from what he told me about it, it is 
perfectly up to date and answers all our present-day needs. So far, 
it has been impossible, unfortunately, to find a publisher willing to 
bring it out. 

October 1944 found Blaton helping to organize the Marie Curie 
Sklodowska University in Lublin. At the outset he aecepted the chair 
of General Physies there and in 1945 that of Theoretical Physies. 
From the Autumn of 1946 he was professor of Theoretical Mechanies 
at the Jagellonian University in Cracow. 

In May 1947, at the first post-war symposium on physies in 
Warsaw, Blaton spoke on the collision of elementary particles in 
relativistie mechanics. He further developed this work during his 
visit to Copenhagen in the Autumn of the same year and the first 
few months of the next. It appeared posthumously in the publications 
of the Danish Academy. 

In Copenhagen, Blaton also devoted much time and thought to 
the problem of the forces causing 4-mesons to enter into the compo- 
sition of a-mesons. He surmised that these forces are caused by the 
electron-neutrino field. One week before his death he lectured on 
these deliberations at the second Warsaw symposium in May 1948, 
emphasizing, however, that his idea was open to criticism. Despite 
the fact that he did not publish his considerations we find a mention 
of them in a note of Professor O. Klein. 

Thus it was that Jan Blaton met his death at the very moment 
when, following the forced inactivity of the war years, he was entering 
a period of increased scientific activity. He recommenced his seien- 
tifie work with wonderful energy. Pitiless fate, however, plucked him 
from our midst at the very moment when he was entering upon the 
path which would possibly have led him to his life's greatest scientific 
achievement. Among Polish physicists his memory will remain alive 
for ever, just as I shall keep for ever in my soul the memory of one 
who was my good friend as well as a pupil dear to my heart. 

In 1949 Jan Blaton was posthumously awarded a National Scien- 
tific Prize. He was being the first Polish physicist to receive it. 


Papers published by the late Jan Blaton 


— 


| (1) Über die Intensitäten der Multipollinien in der Balmerserie, Z. Physa 61, 
263 (1930). 


(1a) O natezeniu linij multipolowych w serii Balmera, Spraw. i Prace Polsk. 
Tow. Fiz., 5, 17 (1930). 


(2) Gibt es eine Doppelstreung von Lichtquanten? Z. Phys., 69, 835 (1931). 
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(3) Uber die Richtungsquantelung der Atome durch eine Lichtwelle, Bull. 
Int. Acad. Polonaise Sci. Lett., Série A, 599 (1931). 


(4) Uber die Dispersion des Lichtes in der Umgebung von Quadrupollinien, 
Z. Phys., 74, 418 (1932); 82, 684 (1935). 

(5) Die Quadrupolstrahlung (with A. Rubinowicz), Erg. exakt. Naturwiss., 
11, 176 (1932). 

(6) The Nature of the Forbidden Lines in the PbI Spectrum (with H. Nie- 
wodniczanski), Phys. Rev., 45, 64 (1934). 

(7) O natezeniach linij dipolowych magnetycznych. Wilno 1934. 


(7a) Uber die Intensitäten magnetischer Dipollinien, Z. Phys., 89, 155 
(1934). 


(8) Quaternionen, Semivektoren und Spinoren, Z. Phys., 95, 337 (1935). 

(9) Versuch einer Anwendung des Fermat’schen Prinzips auf geophysikalische 
Wellenprobleme, Biul. Tow. Geofiz. w Warszawie, nr 14 (1937). 

(10) Zur Theorie der Multipolstrahlung, Acta Phys., Polonica, 6, 256 (1937). 

(11) Zur Kinematik und Dynamik nichtstationärer Luftstrómungen, Biul. 
Tow. Geofiz. w Warszawie, nr 15 (1938). 

(12) On a Geometrical Interpretation of Energy and Momentum Conservation 
in Atomic Collisions and Disintegration Processes, K. Danske Vidensk. Selsk. 
Mat.-fys. Medd., 24, No 20 (1950). 
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ON THE THEORY OF SPONTANEOUS AND „FORCED“ 
K-CAPTURE 


By Marian GÜNTHER, Institute of Theoretical Mechanics, Univer- 
sity of Warsaw, Warsaw 


(received November 27, 1948) 


The main purpose of this paper is to calculate the probability of electron 
K-capture on the basis of the theory of interaction of nucleons with electron-neutrino 
fields (1), in such à way as to avoid as far as possible all special assumptions 
regarding the magnitude of the nucleon radius and the behaviour of the electron 
wave-function inside it. Our final formulas contain only the range of the nuclear 
forees, i. e. the meson mass. 

We shall establish also the dependence of the probability of a K-capture 
proeess on the atomie number Z of the element which is thus transformed into 
another one with the atomic number Z —1. Finally, we shall be led to draw some 
consequences regarding the hypothetical „forced“ K-capture probability when 
irradiating a suitable material by antineutrinos — this being a consequence of 
the neutrino-hole theory. 


Hamiltonian of the problem 
In order to introduce consistent notations for all the five fields 
(electron, neutrino, meson, proton, and neutron fields) it is necessary 
to give first a short survey of their basie equations. 


(a) Electron field 


The differential Hamiltonian of the electron field is 


> > > > 
Ha = (eo + a(ep—eA)-- me jy, -(1) 


where e is the (negative) electron charge. The spinor UN of y 
obey the following anticommutation relations! 


[vo (2), vo ()] , = RR): | (2) 


Introducing the orthonormal set of functions (c-numbers) un describing 
the different states of the electrons and being the eigenfunctions of 


1 As argument of a function x stands for 4,25, 23. 
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the time-free wave equation, according to the formulas 


>> > 
EO u — teo + a(cp—eA)+ Bmw, (3) 
f a umdr= bam; = 


we can expand the y function (q-number) in the series 
= Sanu , (5) 


with the a’s obeying the commutation relations (2) 


[goa 0 (6) 
The total integral Hamiltonian of the electron field becomes thus 
An Hand > dea, — NN HEN ara, (7) 


We must further remember that the application of the a, or y oper- 
ations to the state function means annihilation and the application 
of d, or y — creation of electrons. 

` The operators N describe the number of electrons occupying 
the n-th state with the eigenvalues 0 and 1. 

When dealing with hole-theory it is also convenient to introduce 

the operators 
anīūnr=NP=1—N® for E9<O, ! (9) 


i.e. for the states of negative energy accounting for the number of 
positrons, the energy of these positrons being 


EP=-EN. "etu 


(b) Neutrino field 


All the neutrino field quantities obey the same equations as the 
quantities of the electron field, the only difference being that we must 
put m=0 and e— 0. 

We shall use the following notations for the neutrino field quan- 
tities 

v for the neutrino field, instead of wv for the electron field, 
b, 33 ” 


29 33 » oF An 25 » 55 33 


7 de ) ” »» bj ,5 2) 3) Un 3” ” ” >>) 
n (e) 

Js ” 33 bb) 3) 322. 3” En 39 33 39 29 
Eo Eo 

( ) 55 ” 3) ” ” ” n 3) 3) 3) 3 
n (e) 

N US 393 ,» 39 35 » » NY » LE] 33 39 
NO NO 

v ” $5 39 23 2 33 n 35 3? 33 39 
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Thus, instead of (7) we get 


Hio — f Hatedt= >b,b, B= S NO EQ. (11) 


v 


We shall deal here consistently with the convention that in processes 
like the K-capture electron and antineutrino or positron and neutrino 
are always created and annihilated together. Therefore, only py- or yy- 
like products will appear in the interaction term. 


(c) Proton and neutron fields 


All that has been said about the eléctron field is also valid, in 
principle at least, for the proton and neutron fields, only the wave 
equations Ea to (3)) can be more caret owing to the 
forces of non electromagnetic origin. 

We introduce the following notations 


¥ for proton field, Ø for neutron field, instead of y for electron field, 


B, » 55 33 AN 35: 33 33 33 33 An 3 » » 
Vu » » » Um » » » ? » Un » » » 
N » » » APA » » » » Jes » » 
EY $2 » » ra » » » » » Tuo » » 35 


Then we have, instead of (3), 


for protons Hy= [yir XB.BLES — NER (12) 
u 
W) 


and for neutrons H nour = ji Host dr = Am Rune ANTES. 


(d) Meson field 


The meson field we shall deal with here is a certain combination 
of the charged vector and pseudosclalar fields. We shall start first 
with the Meller-Rosenfeld combination, enabling us to perform the 
calculations for the vector and the pseudoscalar parts of the field until 
a certain point in a uniform way, owing to the formal invariance of 
the equations against five dimensional rotations. To get the other 
combinations after the vector pseudoscalar parts have been separated 
means simply to ascribe different coupling constants to their inter- 
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action terms with the other fields. The equations of the Meller-Rosenfeld 
field are (4) 


Y 9%, 9x, qj q 
Cus = Or, 9% d- 82»; Duy =p, (14) 
4 El 
EICH 2x, 
Dx ug. Ma, (n9 0,1,2,9,4) — (15) N 3-9 (16) 


v—0 v=0 


We can get the equations (14) and (15) together with their conjugates 
from the following differential Lagrangian 


4 4 4 á 4 4 
I=—4)' 2, Fv G ww— D 2 tutu 2; xa My td) xs Mi. (17) 
A—0 v—0 u=0 u=0 u=0 u=0 
To have the terms needed in our problem we must put (5 ) 

M,—g y! yop + uP yor, 
i 2 1 
Sp =F VOL p+ He r9») W. (u,9=0,1,2,3,4) 09 


the g's being the coupling constants. Quantities with the »] -sign are 

defined as follows: 
y'=ipy®=iaß (the same for 9,6, V), 

Gi, Gun Durs nS N=, Ti Se (u,v— 0,1,2,3), (19) 

Gu=—Ou, er —Sya, Mi=—M,, = KA. (20) 

The well-known relations between the y's and the Dirac a and p opera- 

tors are 
y= — ifa, (v= 1 12; 3), OL P, y 9 — y 0,0, ,4) — 1040305. (21) 


In order to split the field into its vector and pseudoscalar parts it is 
sufficient to observe that the fifth coordinate & was introduced from 
a purely formal point of view, i. e. all field quantities are independent 


of it. In this way we get instead of (14), (15) and (16) two independent 
Systems of equations: 


* odd, Q3) 
= v u 
GM ne See (n,v=1,2,3,4) 


4 
9G, 9 
and IX + uM M, | Go, — on +8» (v=1,2,3,4), (23) 


vl 


—— n 
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the eqs. (22) describing the vector (3), and (23) the pseudoscalar part 
of the field. According to the remark made before it is sufficient to 
change the g's in (18) into So and M, to get other combinations of 
these field than the Meller-Rosenfeld. mixture. 

Introducing the canonical variables z, by means of the Lagran- 
gian (17) 


eL diss. + OL 1 
S 2i us ic Cs a a (u=0,1,2,3,4) (24) 
from which it follows that 
a,=n=0, at= n, (v=0,1,2,3) (25) 


and eliminating y, and zł as having vanishing dynamical conjugates 
with the help of the equations (15) and (24), as follows: 


QUE 3 
1 ie yO} I © NELA 
Lig caca eo hae > Ma en ony (26) 
v=1 y—1 
we get the differential Hamiltonian of the meson field and phe inter- 
action in the form 
Fr HO CLHOS-H- HG, (27) 


These four terms correspond to the vector and pseudoscalar parts 
of the field and to the interaction of the vector and pseudoscalar parts 
with other fields NE and are 


(0) E >>. æ > 2 >> 
HO, = rot 7 x: rot jns a ) + Qi a8" n- div zr gàyx, (28) 
Ho = grad y Xo grad phone Wh XoXo» (29) 
ne rot z +8, rob y+ 5, 3-5 <M, divi S dva ; M,M, (30) 
>> >> 
—y7,M — PM Mia 8. Te PS “09 

> » 

Hy, — S.o grad Xo+ | grad yo+ 3.8, (31) 


—ZXo Mo — xo Mo — i62, Int iem S40- | 
Ordinary threedimensional notations have been used in these formulas, 
the threedimensional vectors occuring there, having the components: 
> j E > 3 
X= (Lay Hoy Xs) T= (087.75) M=(M,M,M;) 


S = (S53, S31, 813), S4. (Sai, 845, 845), N.o = (819, S20) Sao): 
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Finally, to quantize the Hamiltonian according to Bose statistics, 
we postulate the usual commutation relations: 


[n(«), xu(&')]- = [n(v), xu(&)]- — —tty Op» 0( —2") (33) 


all other pairs of meson field quantities commuting. 


(d) Total Hamiltonian 


In order to get the total Hamiltonian of all the fields, we must 
add the integral hamiltonians of these tields and those of the inter- 
action terms (space integrals of eqs. (30) and (31)) and — strictly 
speaking — the terms due to electromagnetic interactions also (the 
one, at least, which corresponds to the interaction between protons 
and electrons). But, as we shall see later, we need only know the form 
of the zero Hamiltonian for the meson field and the „meson“ inter- 
action terms given by the equations (28), (29), (30) and (31). 


Canonical transformation of the Hamiltonian 


The „source“ quantities in the interaction terms (30) and (31) 
of the meson field Hamiltonian can be divided into two groups with 
respect to their magnitude in the case of slow ne u: 


(electrons, protons, etc.). Namely, there are quantities as M, M, 5 Sa. N 
the expectation values of which vanish for slow particles -— vanishing 


velocity, and others as M,, 8, Sy with non-vanishing expectation 
values. 

As we shall limit ourselves to nonrelativistie velocities for protons 
and neutrons contained in the nucleus (and even for the electrons 
in the K- mer, we can eliminate, in this approximation, the meson 


quantities a Xo: Fe To from the interaction terms by a procedure equi- 
valent to the introduction of the static Yukawa potential instead of 
them. We can do this by means of a suitable canonical transformation 
of the Hamiltonian. We proceed as follows: first of all, we divide the 
sum of Hiet and Hj, into three parts: 


the ,big" one 


do EEE S 

H®=8 rot y4- rot 5 —— M, divz+— M, diva 
M u 

E (34) 

> ES <= | 

+ S.o grad yot S.o, grad Zo, 
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the „small“ one 


2 Se re n r: 
and finally 
DEN TN >> 


the last one not containing meson field quantities at all. We write 
down the transformation matrix in the form eiX and have for every 
operator Q | 

Q= e'KQe—ik, (37) 


where Y means the operator in the „new“ coordinates. We have thus 


P ma eA (BRIF HO F Hoa Hye) = 
O E OHHO R0. H+H”) 
X(1—iK —4K?--..)— . (38) 
=H wat Hye + iK, Hwat Hi] + H+ ip, HOL 
—4[K,[K, HR+ Hio] J_+H"+... 
If we choose now K to satisfy 


[K,HR tH] =H”, (39) 
we get 
Ana = HOS HO SLE BU] H". (40) 
Now, introducing the Yukawa potential function (7) 
s 
Vina 7 (41) 
which satisfies the „static“ meson field equation 
(A—,2MV(r)- 0 for r0 (42) 


. and remembering the formula 
f V) (4—22s()dv— f V((A—2)o—9(4—,8) V(r)dr——9(0) (43) 


valid for mb functions p, we get for K satisfying (39) the expres- 
sion 


=; Je * Mq) "a xia) ~ 5 Male) div’ a") 
+ > > > 
— Six) rot’ zz) — S(x) rot x(x) (44) 
5. (X) grad’ zu) 8.08) grad’ z,(x))dvd:'. 
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To verify it, we must use the commutation relations (33) and the 
formula (43). In order to show how the calculation runs, we shall cal- 
culate for instance the following. commutator 


|f f TVer) av ple" ded’, n 
— el f fer: )V(xz—2a') div’ M" ‘\dtdt' fa n(x"), aar yde”! (45) 
+5] ffm (x) V (a — x’) div’ M )dvd«' f aiv Mar) div aarde”) 


Taking into account the commutation relations (33) and performing 
some partial integrations, we have 


[/ f a4) V(v—a') div’ aydir, [Reale ae] = 


E SEHR NEAR Völr—a”)drdrde' = (46) 


— ih f f Mia) V(z—a') div’ ae 2')dvdv', 
similarly 


[| [ Male) V(@—2’) div’ a’ )drdr’, f div’ Pa") div” x(0"")dz” | = 


= > (47) 
— —ih sf f V(z— x’) A'{ M (æ) div’ a(x’) }drdr', 
and finally, with the help of (43), 
ES AA 
| 3 M.(2)V(@—a') div’ x(e’)drar', E 
Ee TAN 
= — y?) {M(«) div’ z(z )jydxdc' (48) 


_ the T 
i "E a')dr, 


where the Yukawa potential has disappeared. Calculations with 
other commutators are similar. While performing the partial integra- 
tions mentioned above, we must pay special attention never to introduce 
second derivatives of V into volume integrals or et derivatives 
into surface integrals. 


The only term we have to calculate in our problem is now 
H—5[K,H9] 4 H", (49) 


this being valid also if we perform the transformation taking. into 
account other parts of the total Hamiltonian corresponding to other 
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fields. Indeed, we can easily verify that the commutators of K with 
the zero Hamiltonians of those fields (electron, neutrino, proton and 
neutron fields) and the electrostatical interaction terms between them, 
have their zero matrix element corresponding to the K-capture process. 
The other terms, for which this element does not vanish, are either 
small for non-relativistic velocities of the particles, or can be neglected 
because of their being proportional to higher powers of the coupling 
constants. 

After somewhat lengthy calculations, but proceeding along the 
same lines as in the example given above, we get 


uc ao f (Ve. PARERE Mala) — Mile) Ha’) 


ig EIS, ot” nd o (^ m JSt ) iSt Str) dca 


> > (50) 
-2 f S(x) div’ S(@’) V'Y(z—2a') drdt’ 
+ - 
puff Sol) (V'V(z—2') x rot’ S.q(x')) dvdc' + (...). 
Finally, owing to the eqs. (36), (43) and (49), we get 
H= (fr 0a) (Mo) Mila’) + 42 3. 73 m w')} drat’ 
ae x) div’ Sa) VV (2 —a"') drdr (51) 


TUNIS S.) (V. V(a—a"') x rot^ S, (2) drdr’+(...), 


where the terms marked (...) appear due to non-commutativity of the 
„source“ quantities (as a consequence of the anticommutation rules 
of the form (2)) but have no matrix elements corresponding to our 
problem. It is convenient to express the quantities M,, 5, S. with 
the help of Dirac’s o and c matrices, getting thus 


de E Ed UP) Jo — ay 
M ,— gi yo -r gi T, jane OP rei ER ; 
52 
ow f: op — D oor } 
u u 
where * t: 


We must bear in mind that we may also change the coupling constants 
Jo) 92 ior S.o in order to get a more general field than the Meller-Ro- 


senfeld mixture. 
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The transition matrix element 


The matrix element H corresponding to the K-capture process 
is easily seen to be 


H= f V(z—2a") (g,g, U (x) V (x) Jula") 
2M E = > > 
+ 9,9,(U(a)oV (x), ö&’)ou(®’) + U (x) ogoV (x), d(x )ogou(")) jdvdv' 
Lif f vp RM ln c S 
tog f ST) div’ (5(27) eo (o) V Va) dede 


1 ı TT >” PLE y t= xd "AM , 
—3J f 2s U(x)oV (x), (V'V(z—2') xrot' [v(z') au(z')]) dvdz'. 


The functions u and v describe the neutron and proton states, respec- 
tively, the latter superseding the former in our process, u is the electron 
function in the A-shell, v that of the neutrino plane wave. In our 
further ealeulations we shall neglect the last two terms in (54). We 
ean do this — as shown in the appendix — by assuming the ratio of 
the nuclear radius a (being defined, e. g., as the mean extension of 
the u and v functions) and the range of nuclear forces 1/u to be great 
in comparison with the fine structure constant multiplied by the atomic 
number: 

Za (ay. (55) 


It might also be shown that when dealing with these terms one would 
be forced to make special assumptions about the behaviour of the 
eleetron function inside of the nucleus which we shall avoid. 

On the other hand, in order to be able to perform the further 
calculation without special knowledge about the behaviour of the 4 
and v functions, we must assume the range of the meson forces 1/, to 
be great in comparison to the nuclear radius a, so that we must complete 
the inequality (55) to 

Za««an««1. (56) 


If the last condition can be fulfilled we may write 


mu un 


Hy [V9){9,,(Döz)u(z)+9,93((0,0), B(a)0,0%(@)+(0),0(@)ou(a))}dr, (57) 


> : 
where the radius vector r is to be taken from the electric centre of 
the nucleus, and 


Q)—f'Ut)Y(G)dr, (s) fU(z)sV() de, (5) f Utz)osV (six. (58) 


LE zz 
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For the more general meson field than the Moller-Rosenfeld field we 
should have 


—_— M 


> 


a> EN > > > 
H k— [V(r {gg (Dwu) 9o95( 050), BE oult)+gg0), Loule) jdr. (57 


Nm! 


It is understood also that the approximation (57) or (57) is valid if 
the neutrino wave length is great in comparison with the range of the 
meson forces l/u, i. e. if the neutrino energy is small in comparison 
with the meson rest energy. 

We must insert now for w the solutions of the Dirae equations 
for the K-shell, for v the neutrino plane wave, and then average the 
square of the absolute value of H with respect to: 

1) two different orientations of the spins of the electrons in 
the K-shell, 

2) two different orientations of the neutrino spin with respect 
to its motion, 

3) all possible directions of the rejeeted neutrino, 
as the transition probability depends only on this mean value. 

The Dirac equation of an electron in the central electrostatie 
field of the nucleus of charge —Ze is 


(47 + ae aa) ihoV— jme] u= l. (59) 
If we use the customary representation for the a; and f matrices: 
0001 03030 — 
wep BOO ARO a 3140. 50502 10 
e c loo prot M as no. 7 anos 
1000 1-0 0 0 (60) 
0010 100 0 
OOO hei: o wo <0 
SERM (sO Oe 00 TOT 
017050 CO Or) 
we get the following form of the two solutions for the A-electrons 
1 ya) = 
(D. n —ij, 9 
May; — en 2 
f u® — f'(r) 
uf? — 0, Vix 
1 or T 4 (61) 
u$) = M co80, j m iret sin dep, 
1 
ul) = )-—— sin de, u® — —f(r cos 9. 
fes > 8 —! fh K ! Vas 
2 
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These solutions belong to the common lowest energy value - 


Ex= me V1 zæ, a= ehe (62) 


and correspond to the two possible orientations of the spin of the 
electron, so that if we introduce the angular momentum operator 


> > > à 
m= —ih(r x V)4- tho, (63) 

we will have for its 2-component: | 
mW = $U, m; U—sj,= —4u_y,. (64) 


The functions f(r) and f'(r) satisfy the equations 
(++ meii =o, 
(65) ` 

era) 


and are 


jm pe Vrai re zu 
2 RD 
f=clı+fi- z& CE zt 


with the normalization factor 


(66) - 


Zame VEZEH TA mo Di ge Pe. ae Sg f 
c=(2 = ) Te, ow | (67) | 


The neutrino function v satisfies the equation | 


eS Pun T2 ) 
Y C OEN ec VE 


"bes Pi 
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neutrino energy, and two eigenvalues equal to —1, belonging to negative 
energy. We have thus 


> > 
Ey=—EOQ)=kelk|, EI—=-—EM)=helk|, (71) 
and 
>> >> 
eg D = (a, ko) er, ro) (a, ke, (72) 


where E and E mean the neutrino and antineutrino energies respec- 

tively, the subscripts „+“ and „—“ being easily understood to refer 
to positive and negative energies. It is useful for the sequel to write 
the eqs. (71) in the form 


>> >> 
HAF (a,ky)} = 0, Flle, ko) jeD «D. (72°) 


Inserting (61) and (69) into (54), we get 


Key Ilesha tOO) VOF Fude (T3) 
If, as mentioned above, 


=> 
H»»|k|, i.e. hye>> EO, (74) 
we can write 


DET 1 e—ur—ik a 1 eur 
3 —ik cea. - ee nn x IL / A 7E 
ar i k x ydt= zS x udr zS : udt (75) 


If we further assume, as may safely be done, that the meson mass 
.is great in comparison with the mass of the electron 


hue y» me, (76) 


we get finally after elementary calculations with the help of eqs. (61), 
(66), (67) 


ur 
1 | t udt=Ny where 
4 


77) 
SS v-ze+r IN —-Ze+1) 
Za? | 2Za—~ 

Mig Ten ze s i Verefi-ze+1) Z*a? - 1) 


and where (cf. (60)) 


/ 0 1 
Bu 0 
salo) xo] 2 (09) 
0 0 


9* 
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are the unit „eigenspinors“ of the f-operator, belonging to the eigen 
value 1. If we introduce formally the two other „eigenspinors“ of this 


operator, namely 
0 
0 
J (78) 
0 


U1F 0) =0, F(1t 99) C9 — 7V, (79) 


two equations which are similar to (72’). 
Instead of (73), we get 


we shall have 


W bai f R Saat Oe s 
Hx N RR, eso t (0), 03277 (30) 


Now, we are able to perform the above mentioned averaging of the 
square of the absolute value of the ie (80). We denote by 


a 


Q= js (en + 99: (059), PEN any (81) 


the most important operator in (80), which is of the type of Dirae's 
a and p operators, and average with respect to the neutrino spin orien- 
tations first. This average value is (owing to (72’)) 


4S (4) ELI 18,09 (xC0Qe- Q4 C0) 
>> 
—48.(?Q9Q40 4( (ako) )e:eQ 4C?) — $4609 $(1 ok (a9) QP 


S,(+) denoting the summation over positive or negative energies only, 
while S, over both of thern (we reserve the summation over «C2 for the 
case of „forced“ K-capture). Now we can easily average over all possible 


(82) 


>> 
k-direetions, by simply dropping out (a,k,). 
Denoting the angular mean value by <...... >, we get thus 


GE S,0IH |» — 130299 409. (83) 


Finally, we perform the summation over both possible orientations 
of the electron spin in the K-orbit, and get in a similar way as before 
thanks to (79) 


Hr? — 8,70 (4 8,Co|H «5 — Tr( QQ(4- 03)) (84) 


but without the factor 4 before Sy, owing to the fact that there are 
two electrons on the K-orbit. Inserting the expression (81) and cal- 
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eulating according to well-known rules the „trace“ of the expression 
(QQ(14-0,)), we get 


9 al NOT ONG P WS MES ES > zs OE 2|/1—ze24-1 
|H}? = Saul (19)? (1 + (9595? |(050) + (o)|?} [22a i) 


LL 
, Ii —ze rg ze 2) 
ar(2J1—Z263-- 1) - 


(85) 


To deal with more general combinations of the vector and pseudo- 
scalar meson fields we must replace the bracket expression by 


“unn Ne 
> 


(9.9)2| CDI 4- 19.936059) + 8595(0)) (85) 


Transition probabilities 


In order to get a suitable survey of the possible phenomena it 
is convenient to deal shortly with the perturbation method for the 
case of resonance absorption. The basic equations in our case are (8) 


ine rd II>*<I|H|P>exp [e (Hr—Hn), 


iS = |PICPIH ID exp f 119] + 
it (86) 
+ 3 [F o*<F|H|I> exp I (HH), 
FI 


«s [Fit DHIE exp {Ee Hy, 


where we have used Dirac’s bracket notation (thus, 4...» does not 
mean angular average, as in the preceding section), the suffixes P, I, Fr 
refer to primary, intermediate, and final states respectively, the sum 2 

FT 


means summation over all Fy with fixed I. 

The |P>* state corresponds in our case to an atom of atomic 
number Z with 2 electrons in the A-shell, and to no neutrino in the 
positive energy range. We may postulate further in this |P>* state 
the existence of antineutrino radiation in space, corresponding in the 
neutrino-hole formalism to some unoccupied states of the neutrino 
negative energy range. In this way we shall be led to the aforesaid 
forced" K-capture in case of a suitable (negative) energy difference 
between the initial and the final energy of the atomic system. 

The | I>* states correspond to the atomic system after the K-capture | 
process and one additional neutrino more than in the |P5* state. We 
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call this state „intermediate“ only from à somewhat formal point 
of view, owing to a possible „reabsorption“ of the neutrino (,,re- 
emission“ of an antineutrino). 

The |F7>* states correspond to this „reabsorption“, with the 
atomic system back again in the primary configuration. It is neces- 
sary to ascribe to these states two indices, # and I, the first indicat- 
ing the neutrino which has been absorbed, the second the emitted one. 

The approximate solutions of the equations (86) for |P>* and 
|I1>* are 


o 
[P> = 3! 
i y Prom 21 (87) 
e n € —e 
|I>*=<P|H|I> —— 
Hr—Hı+ ihi 
where | 
27 i 
y =F (KHF oP Lr) pne np; (88) 1 
í 
and where o(Hr,) is the final state density per energy interval, while 


EU Xu uni. DIT e 
r= PA ce 


ES 


which gives: ! ES 


(a) PEE PHD PHD} apap» en) 


m 


"NS 
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We ean now see that: 

(1) In the positive energy range of the neutrinos the density of 
the intermediate |I5* states is the maximum possible one, owing to 
Pauli's exclusion principle: 

OE RISO ha ee 

e(H ) = e,= — cem N (91') 
where Q is the periodicity volume, whereas the final-state density 
e(Hr,) vanishes: 

e(H pj) — 0, (91) 
corresponding to the fact that only one neutrino can be reabsorbed, 
namely this one which has been emitted first, and this leads exactly 
to the restitution of the primary state |P>* only. 

(2) In the negative energy range of the neutrinos the interme- 
diatestate density o(H7) is determined by the intensity of the assumed 
antineutrino radiation, so that if we designate by I(E(?) the number 
of the antineutrinos of energy between E@ and E@+dE® passing 
in unit time perpendicularly through unit area, we will have 


er) SH). (92) 


The final state density o(Hr,)is now at its maximum, because of the 
presence of neutrinos occupying the negative states of energy. Thus 

e(H pj) — Pp: (927) 
In order to know with which of these two cases we have to do, we 
must realize that 


Hp—H)=(E®+EQ)—(BV+E™) or —(EO-EES--EO)— E, (93) 
where EY is given by (62). The condition Hp—H,=0 leads thus to 
po — RO (poa EM) + E= AQ. (94) 


AQ meaning the mass difference measured in energy units between the 
atomic weights of the Z and the Z —1 isobars (nucleus 4- electron shell). 

In the following formulas we may insert |H x|? from the preceding 
section instead of |«I|H |F pl? and |CP|H I»? and introduce an important 
new quantity, namely the transition probability corresponding to 
maximum density of states: 

Qn Eo PN ales I, 
Dp— |AyPep= n ey (II) COP gles o)P + 995 (0)73 
<f Dh mue: rV Z3 +1) VIE +2) 
hu 2r(2]/1—Z*a: +1) 


wv 


(95) 
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Thus, we are led to the following conclusions: 

(1) „Spontaneous“ K-capture. In this case we have 4Q>0, 
i. e. the energy of the atom with atomie number Z is greater than that 
of its isobar of atomie number Z—1 in which it is transmuted. It is 
quite obvious that it is only this that we observe in all the natural 
K-emitters. As 4Q=E™, we have to do with real emission of the 
(positive energy) neutrino and may write 


elle, ee (96) 
The mean lifetime of such a K-emitter is given by 
in pe (97 


(2) „Forced“ K-capture. The energy condition would now demand 
AQ«0, i.e. the isobar Z—1 to possess greater energy than that of 
atomic number Z, and the gap to be filled up by the absorbed anti- 
neutrino. We have now 


y=l, (98) 


but for /' we must still distinguish between the two cases of „wide“ 
and „narrow“. antineutrino lines: (a) „wide“ line 


2 I(E)O IE 
r=% Hf ae E (99,1 (1) gedit e ale) 


(99) 


? 


x (oz a rn Zet r(J1—z268-4- 1) r(J1—22a2 + 2) 
hu 2r(2]1—Z2a24- 1) 


where I(E) is the number of antineutrinos of energy in the interval 
(E,,E--dE) passing in unit of time unit of area perpendicular to the 
direction of their motion. We get as before for the mean lifetime 


QE TTA: (97°) 


(b) „narrow“ line. According to (8 9”), (95) and (98) the transition 
probability per unit time is 


kr): 
20 2 

BN v4 AO 
Yan: re B (=F) 


R(D)— (100) 


Taking into account that the total number of antineutrinos passing 
in unit time unit perpendicular area is Oc/Q, we get thanks to (91) 
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and (93) the following expression for the total cross-section for a single 


antineutrino 
h PA 
A? D 


OR— OLI 1 9 ? (101) 
der SEPA à SNe 
el 


where 4 is the wave-length of the antineutrino, and Ip is still given 
by (95). 


Discussion of the orders of magnitude (9). 


We shall compare first the order of magnitude of the gg'-coeffi- 
cients, which we get by inserting the known experimental data for several 
K-emitters into the eqs. (95) and (97), with those given by Yukawa. 


We assume (1), Tus (e to be of order one, but this order can be 
smaller when the transition is forbidden by nuclear selection rules, 
so that we may compare the short life A-emitters only and establish 
the upper limit for the gg'-estimations, as they will be underestimated 
in general. We must further pay special attention to those short life 
K-emitters which possess small A's, because, owing to (95) and (97), 
the mean lifetime should then have the tendency to decrease with 
increasing Z. While performing these estimations, we assume that 
the meson mass is 200 times that of the electron. 

As a first example we choose the A-reaction 


1 Be 3Li. 
The half lifetime 7 is in this case 53 days, so that the energy difference 
AQ between {Be and $Li, as calculated from their mass defects, is 
to = Un2 = 6,611 08s. 


and we have for the bracket expression in (95), with sufficient approxi- 
mation, 


— 


XA ota Qa, Sukie? 1 (200 
U?) P+ 1929s (050) + 9595 (0j © (AQ)? caes 


= igmbsensedozph 


8 
x 137] 


The g and g’ coefficients estimated by Yukawa are given approxi- 
mately by , - 
(g/bm? 1 


m 10! SAU, 


so that 
(gg P= (2x)* x 5x 10, 
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which gives the same order of magnitude as our estimation. Because 
AQ has not been measured for other K-emitters, we will only give 
one more example, namely that of the reaction 


9 Ga 25 Zn. 


Sp, 
In this case we can be sure that A<<2me? as there is no f radioactivity 


accompanying the K- aue The half lifetime v of $ Ga is the shortest 
one among all „pure“ K-emitters and equals 15 min. Thus 


los) in2-—— 129 10* 8; 


If we earry out the estimation in the same way as before, putting 
simply 4Q9=2me, we get 


je M 


(gg P CDI + Ios go 050) + (9595) (DI «(20 x 27 x 10-5, 


which gives the same order of magnitude as for 1 Be. 

The probability of „forced“ K-capture — if such a K-capture 
exists — would be in general many orders of magnitude smaller than 
the probability of „spontaneous“ A-capture. 

If the formula (99) applies — this would be the case, e. g., when 
irradiating a suitable element by the broad spectral line of antineu- 
trinos produced together whith electrons by a f emitter — the ratio 
of probabilities of „forced“ and „spontaneous“ K-captures- will be 
of the order of magnitude of the ratio between the radiation density 
owing to the presence of these emitters to the maximum state density 
provided by Pauli's principle. We shall estimate the order of magnitude 
of this ratio in the centre of a sphere of radius R containing the uni- 
formly distributed emitter, which we can choose to be radium C pos- 
sessing the half lifetime 


TR= 19,5 min., tg — v g[ln 2 
and the maximum f-ray energy 
Emax— —3,15 MeV. 


After elementary caleulations we get for the intensity in the sphere 


3L 
I= 2131088. 


1672 tr Emax EB R (R in em), 
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where L is the Loschmidt number and » the number of gram molecules 
of Ra contained in the sphere. If we calculate now the maximum 
energy density (91), assuming for instance E=2me, we get 


Ir/Imax= 6 x 10-23 y[R2, 


where Imax is the radiation density corresponding to the states density On 

If we assume two elements — one of which is a natural K-emitter 
and the other the hypothetical „forced“ K-emitter — to have equal |H g|? 
and equal absolute values |4Q|=2mc and place v’ moles of the latter 
in the centre of our sphere, we get for the number of these „forced“ 
K-processes in unit time 


t being the mean lifetime of the natural K-emitter. 

I wish to express my deepest thanks to Prof. A. Rubinowicz, 
to whom I owe the idea of this paper, for his suggestions and constant 
interest during the work. 


Appendix. 


We have still to discuss the order of magnitude of the last two 
terms in (54) which we have neglected. 

Using the same approximation as in (54) (based on the assump- 
tion (55)), we get for these two terms 


"T vere en > 
520), f ET div (2trjas sur) de 


and (a) 
1 Ax V(r) > E 
— tetto), f x rot {oer)ou(r)} dr 
' respectively, where V,(r)—dV (r)/dr. 
We can easily see that 


> = ah eA 
div (v(r)esou(r)j — v(r)es(o)Vu(r) + u(r) o5(o V) (7) = (b) 


2 EU. tap WAS 
= 10(r) 0,4 Vu(r) — iu(r) osa V(r). 


Thanks to (56) and (69), this is equal to 


sre o s — is ua ona) u(r) 
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which, strictly speaking, would even lead to divergent integrals, owing 
to the term Za/r, if we assumed a too small nuclear radius. In order 
that the last expression be small, for every r, in comparison with 


uö(r)u(r) and like terms (c) 


we must put Za/a<<u, thus giving us the condition (55). 

Performing the integration in (a) shows that when div {öfr)ozou(r)} 
is small in comparison with the (c) — terms the first term in (a) is 
small in comparison with the other terms in (54) and can thus be 
neglected. — i 

In a similar way, we get the same condition (55) ) when the second 
term in (a) has to be neglected. S 
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A NOTE ON THE INVARIANT FORMULATION OF THE 
QUANTUM FIELD THEORY 


By Jerzy RAYSKI, Nicholas Copernieus University, Torun 


(rece ved June 18, 1949) 


The Schröd.nger equation is written in a relativistically invariant form in 
the general ease where the interaction-energy density is not a scalar. The invariance 
may be exhibited without introducing a generalization called ,infinitely-many-time 
formalism". The autonomy of the formulation of Tomonaga, i. e., its independence 
of the Lagrangean and Hamiltonian formalisms, is emphasized. 


The famous relativistieally invariant formulation of the theory 
of quantized fields developed by Tomonaga. (1946) consists of two 
decisive steps: (a) The transition to-the interaction representation 
where the field variables obey interaction-free field equations and 
covariant commutation relations. (b) The rewritting of the Schrö- 
dinger equation in an invariant form: 


OVS] = 


i 
ÒS p 


H(P)PLS]. (1) 
(1) is not yet sufficiently general as it applies only to the case of a scalar 
Hamiltonian density of interaction. In general, the operator H is not 
a scalar but may be represented as the four-four component of 
a tensor I,,. The special case of a scalar H is contained in the more 
general one if Iw is of the form — Hoy. 

An interesting feature of the Tomonaga theory is that it permits 
to forget the Lagrangean and Hamiltonian formalisms: The field 
equations and the commutation relations in the interaction repre- 
sentation are sufficiently simple to be regarded as primary concepts. 
They are understandable by someone who has never heard, e. g., about 
the canonically conjugated momenta. The tensor Iw may be postu- 
lated (similarly as there was always postulated a Lagrangean) and 
we are not obliged to remember that its four-four component is identical 
with (the interaction part of) a Hamiltonian. This statement seems 
to the author not quite trivial since there exist possibilities of theories 
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which are not derivable from a Hamiltonian formalism but which 
might be expressed in terms of the Tomonaga formulation !. 

The replacement of the planes t= const. by more general space- 
like surfaces (i.e. the so called intinitely-many-time formalism) has 
been introduced in order to give the Schródinger equation an invariant 
form. However, it is evident that in the frame of the restricted relativity, 
it is not necessary to use this last generalization in order to write the 
Schródinger equation in am invariant form. In the theory of special 
relativity a space-like plane 2 is also an invariant concept so that 
we do not need to introduce arbitrary surfaces. 

We consider a functional V[X] and define a „normal derivative 
of V[X]in the plane &,“. We take a set of planes 2; parallel to 2, 
and form the differential quotient 


V[24] — PLZ,] 
l ? 


where | means the four-dimensional distance between the hyperplanes 
2; and Zo. If the limes exists for 2,—2, it will be denoted by 
dip ui 
dX, 


and called „the normal derivative in the plane Z,*. 
Given is a tensor I,,, then we may write the equation 


dV[X| : M 
dx =ef daal ot LES, (2) 


X 


a 


where 5 do, denotes the surface integral over the plane 2 and n, is 
. 54 


a unit vector normal to the plane X (n,m,— —1). For definitness, we 
may also assume that n, points to the future, (ny — 0). (2) is an invariant 
equation since it makes use only of geometrical concepts independent 
of any special frame of reference. 

If the system may be expressed in terms of the Hamiltonian 
formalism then Jo, will be identical with the interaction part of a Ha- 
miltonian and (2) will be equivalent to the Schródinger equation in 
the interaction representation. To see this, it is sufficient to write (2) 
in a special Lorentz frame in which the planes t= const. are parallel 


1 As examples of non-Hamiltonian systems we mention the theories where 
the interaction is cutt off (spread out) by means of integral (Peierls & Mc Manus) 
or differential operators (J. Rayski, 1947). 


al 


Invariant Formulation of Quantum Field T heory 31 


to the plane 2,. In this system of reference the functional V/[ X] may 
be replaced by a function of t, while the normal derivative becomes 
the usual derivative with respect to t. Moreover, nuIy,n,—I,— H so 
that we get from (2) the traditional Schrödinger equation 


EDAC | 
i‘ = "= f esi. (2^) 


By comparing (2) with (2') we recognize the double meaning 
of the variable ¢ in the Schrödinger equation. In H(æyzt) it appears 
as a space-time variable, while in Y(t) it is only a substitute for the 
four-dimensional distance between two parallel space-like planes. 

Of course, the restriction to planes may be as well abandoned. 
It is easily seen how (2) should be generalized for arbitrary space-like 
surfaces S: 


SPIS 
re = my(P) L(P)n,(P) PES]. (3) 


In case 1,,— —H 0,, the equation (3) becomes the equation (1) of Tomo- 
naga. Let us denote by C that part of S+ ôSp, which does not over- 
lap S. The points on C will be denoted by Pc. For reasons of consistency 
it is required that (a) C contracts to the point P, (P5— P$—-0) and (b) the 
angle between any normal to the surface C and the normal n, (Po) 
tends to zero: n4(Pc)—nQ4(Pg9) —0 for 05—0. Due to these conditions 
the variations are weak and the functional derivative ö/öSp, refers 
to and only to the properties of the surface S in the vicinity of the 
point P,. 

A formulation equivalent to our (3) has been given by Matthews 
(1949) by means of a generalization of the Hamiltonian: H —n,I,n,. 

(2) is a special case of (3), namely, we may replace in (3) the 
surface S by a plane & and, by means of a set of successive local varia- 
tions 6S, we may perform a transition to another plane &+d2 parallel 
to X. By summing up the local variations we obtain in limes the equa- 
tion (2) ?. 


i 
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INFLUENCE OF TORSIONAL VIBRATIONS OF LUMINES- 
CENT MOLECULES ON THE FUNDAMENTAL POLARI- 
ZATION OF PHOTOLUMINESCENCE OF SOLUTIONS 


By Alexander JABLONSKI, Physics Department, Nicholas Copernicus 
University, Torun. 
(received June 19, 1949) 


The observed values of fundamental polarization are always smalier tuau 
those to be expected on the ground of theoretical considerations. This fact can be 
explained (at least partially) by the influence of torsional vibrations of fluorescent 
molecules on the rate of polarization of photoluminescence. Some conclusions can 
be drawn from the difference between theoretical and experimental values of the 
polarization. This note contains some provisional results of theoretical investigation. 
A fuller report is intended to be published shortly. 


According to Pauling (1930) and Stern (1931) the molecules (and 
parts of molecules) in many crystals can rotate if the temperature 
of the crystal is sufficiently high. At low temperatures the molecules 
carry out small oscillations (torsional vibrations) about their equi- 
librium orientations. The energy levels and the corresponding eigen- 
functions approach in the case of diatomie molecules to those of a two 
dimensional oscillator. The case of polyatomie molecules (three finite 
principal moments of inertia) was studied by Mrs W. Hanus (to be 
published shortly). In the last ease the levels and eigenfunctions of 
the lowest torsional vibration states are very nearly those of a three- 
dimensional harmonie oscillator. Thus molecules in crystals must 
possess an amount of torsional vibration energy even in their lowest 
state („zero point energy“). 

There is no doubt that similar torsional vibrations are carried 
out by luminescent molecules in solid (vitreous) or very viscous solutions. 
These vibrations (and sometimes also some of the normal internal 
vibrations of molecules) must cause a partial depolarization of photo- 
luminescence. Since this cause persists down to the lowest temperatures 
(zero point energy!) it cannot be thoroughly eliminated by choosing 
suitable experimental/conditions 1. Hence, the explanation of the fact 
^ 1 In contradistinction to the two other thus far known causes, i. e. to the 
Brownian rotation of luminescent molecules and the transference of the excitation 
energy from one molecule to another. 


Acta Physica Polonica 3 


34 Aleksander Jablonski 


that the observed values of the fundamental polarization ? are always 
smaller than those to be expected on the ground of theoretical con- 
siderations. 

A theory of the fundamental polarization for molecules at absolute 
rest (not executing torsional vibrations) was given by the present 
writer (1935, 1936). According to it, the fundamental polarization for 
the ease in which the same virtual electronic oscillator is responsible 
for the absorption as well as for the emission of light is given by? 


3 3 
3Xri-(Xrny 


Hi E] 
De a (1) 


XresQuy 


i=1 


where J, denote the relative principal transition probabilities along 
the three mutually perpendicular (principal) axes of the virtual electronic 
oscillator. 

From (1) one obtains P,—1/2 for a linear oscillator (1414-0, 
I,— 103-0), P,=1/7 for a symmetrical flat oscillator (JL, — 54-0, 75— 0) 
and P,=0 for a spherical oscillator (/,—5— praem 

Let us now take into consideration the influence of the torsional 
vibrations of luminescent molecules on the rate of polarization. In this 
note we restrict ourselves to two simplest cases. Let y denote the 
angle of deviation of one of the principal axes from its equilibrium orien- 
tation and let the average value of cos? y be cos? y=1—sin? y=1—e. 
Provided the distribution of y be axially symmetrical, a linear virtual 
oscillator becomes (owing to the torsional vibrations) equivalent 
to an axially symmetrical spatial (three dimensional) oscillator with 
relative transition probabilities F, = 1—e, I,=I,=e/2 (instead 
of J,=1, l,—10,—0; the common factor is ommited as irrelevant). 
Putting these values into (1) we obtain 4 


r 0225312614 


Thus Po<1/2 for e» 0. 


? The limiting value of the polarization of luminescence of an isotropic solution 
observed at right angle to the electric vector of plane polarized primary light when 
depolarizing factors (footnote (-)) are eliminated. 

* This formula is given in a different but equivalent form in papers Jabłoń- 
ski (1935) and (1936). 

* This formula is equivalent to a formula derived by F. Perrin (1929) for 
a linear oscillator carrying out quick irregular oscillations about its mean orientation. 
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For an axially symmetrical flat oscillator the transition proba- 
bilities become (provided the distribution of y for the symmetry axis 
perpendicular to the plane of the oscillator be axially symmetrical) 


Je Sin? y— es, 1,1, Pa (instead of /,=T,=1, T,=0) and hence 


by means of (1) 
92 — 12244 


Pe cp - 
9 32—4e4-28' (3) 


which shows that Po <1/7 for e>0. 

We should like to postpone the publication of more general for- 
mulae to a later paper. For the time being we admit an approximate 
validity of (2) and (3) for particular cases discussed below and draw 
some conclusions from experimental data published by Perrin (1929) 
and by Feofilov (1947). 

Assuming the virtual oscillator to be linear in the case of fluores- 
ceine and resorufine and using the experimental value Po=0,44 of 
Perrin, we calculate e by means of (2) We obtain e= 0,048. Provided 
the whole „polarization deficiency" P,—P$—0,5—0,44— 0,06 is due 
to the torsional vibrations of the above molecules, there results a stan- 
dard deviation of the direction of the virtual oscillator from its equilib- 
rium orientation equal to 10° for fluoresceine and resorufine solutions 
in glycerine at room temperature. 

Using Feofilov's value P; — 1/14 for the fluorescence of benzene 
solution in glycerine and assuming the oscillator in this case to be flat 
and symmetrical (the transition moment is known to be parallel to 
the plane of the benzene molecule for the band under consideration), 
we obtain by means of (3) £— 0,2, which corresponds to a standard 
deviation of the sixfold symmetry axis of the benzene molecule equal 
to 26°. This value seems to be somewhat to large. Possibly there were 
additional depolarising factors in Feofilov's experiments. 

In general the distribution of y depends on the principal moments - 
o inertia, on intermolecular forces, and on temperature. The polari- 
zation deficiency depends apart from the above factors also on the 
anisotropy of the virtual oscillator and on its orientation with respect 
to the principal axes of inertia of the fluorescent molecule. Postponing 
the discussion of the general case to a later paper, we confine ourselves 
to a rough treatement of a simple case of a flat molecule having an at 
least threefold symmetry axis and a transition moment lying in the 
plane of the molecule, i. e. to the fluorescence of the benzene molecule 
in glycerine solution. In this case two of the three principal moments of 


inertia are equal. The third moment of inertia is irrelevant. We assume 
3* 
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also the two relevant torsional eigenfrequencies to be equal. Then 
the angular frequency of this vibration is o=Yr/I, where v is the mo- 
ment of torsion and I the moment of inertia. Assuming the eigen- 
functions to be those of a two-dimensional symmetrical harmonie 
oscillator, we obtain the following formula for the mean square of y 


— 2h 1 1 
Joc zn 
7 ler +a) 


erT — 1 
for kI»»ho 


Iowo Iœ?’ (4) 


Since e— sin? y and sin? ywy? for y<<1, we can easily calculate œw on 
ground of the values of e deduced from experiment by means 
of (3), provided the moment of inertia of the molecule is known. Put- 
ting 1—1-4.10 ?g.cm? and T= 263? K, we obtain ¥ —o[2zc— 26 em-1 
as order of magnitude of the frequency of torsional vibrations about 
axes lying in the plane of the benzene molecule for benzene solution 
in glycerine. Thus the frequency appears to be of the same order of 
magnitude as that deduced from the Raman spectra produced by 
torsional vibrations of different molecules in crystals (Rousset, 1947). 

As can be easily seen e does not vanish when the temperature 
of the luminescent solution tends to 0? K. At 0 K 


TEE 
Io YI 
It seems probable that the-investigation of the fundamental 
polarization and of its dependence on temperature can provide us in 


some cases with informations concerning intermolecular forces. 
A fuller report is intended to be published shortly. 


References 


Feofilov P. P., C. R. Acad. Sei. URSS., 57, 343 (1947). 

Jablonski A., Z. Phys., 96, 236 (1935); Acta Phys. Polonica, IV, 371, 
389 (1935); ibid., V, 271 (1936). 

Pauling L., Phys. Rev., 36, 430 (1930). 

Perrin F., Ann. Phys., Paris, 12, 169 (1929). 


~~ Rousset A., La diffusion de la lumière par les molecules rigides, Paris 
(1947). 


Stern T. E., Proc. Roy. Soc. A., 130, 551 (1931). 


Vol. X (1950) Acta Physica Polonica Fase. 1—2 


THE PHENOMENA OF MOLECULAR ORIENTATION 
IN POLAR LIQUIDS AND THEIR SOLUTIONS. PART I. 
EXTENTION OF ONSAGER'S THEORY 


By A. PIEKARA, Institute of Technology, Gdansk. 


(received July 8, 1949) 


$1. The assumptions of Onsager's theory are summarized; the view that 
all molecules surrounding a given molecule may be considered as a continuous medium 
is rejected as inadmissible. $ 2. Onsager's theory is generalized by treating the 
nearest molecules in the neighbourhood of the given molecule individually, whereas 
all those further away are considered as a continuous medium. The Onsager sphere 
separating the nearest neighbourhood of the molecule from the continuous medium 
has to be increased in size. $ 3. Formula (36) is deduced and by means of this the 


molecular dipole polarization pur of a given polar substance, pure or dissolved, 
- : s , o 4uNu 
may be calculated. The ratio of this quantity to P= Tu: 
duction factor R. The dependence of zur and E on the concentration of nitrobenzene 
dissolved in benzene is shown graphically (fig. 1). It is evident from Onsager's curves 
that for greater concentrations a new kind of coupling appears increasing the polari- 
zation. It consists probably in the coupling of molecules into ensembles greater than 
pairs. It is possible that antiparallel pairs are coupled in parallel (fig. 2c). $ 4. For- 


mulae (43) and (64) for the modified Kerr constant K"" are deduced and applied 


to solutions of substances in non-polar solvents. Hence, the molar constant Kz of 
the dissolved substance may be calculated. The results for nitrobenzene are pre- 
sented graphically (fig. 4). Here also Onsager’s curve in contrast to Lorentz’s curve 
rises for greater concentrations which seems to suggest, as in the case of polarization, 
a multiple coupling. § 5. The increase of the dielectric constant Ae under the influence 


is called the re- 


of an external electric field is expressed by formula (84). The molar constant S" of 
electric saturation is defined in (86). For dipolar substances dissolved in a non-polar 
solvent they may be determined from formula (85). Experimental data for the same 
nitrobenzene solutions are illustrated graphically in fig. 5. § 6. In the Cotton-Mouton 
effect the orientation of molecules is due-to an external magnetic field. Owing to the 
extremely small magnetic permeability of organic liquids the external and the local 
fields are identical. For this reason Onsager’s idea cannot produce any different 
results in this case. The Cotton-Mouton molar constant C" as defined by (92) 
is expressed by (93) and (94). For liquids formula (96), expressing the additivity of 
the constant €", is applied. The molar constant Of for nitrobenzene calculated by 
means of the latter formula is presented graphically in fig. 6. Its increase with in- 
creasing concentration suggests coupling of molecules in aggregates of two or more. 
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8 1. The fields of Lorentz and Onsager. When in a dielectric 
medium there exists an electric field of intensity E, then each molecule 
is under the action of a field F, which is a little stronger than E, and 
is called molecular field. This is due to the fact that the field E is 
ex definitione a field acting inside a long and narrow tubular cavity 
cut in the medium along the lines of forces, while the molecules of 
this medium are subjected to an additional field due to the surrounding 
polarized molecules. 

Lorentz has calculated the molecular field in a well-known way. 
He surrounds the molecule under consideration with a sphere from 
which all other molecules are removed, assuming that this does not 
alter the field F, nor the distribution of force lines inside the sphere. 
The field is given by the equation 


F=E+ 9, (1) 


4nJ|3 (J —the electric moment of 1 c. e. of the medium) is the addi- 
tional field due to the electric charges of the polarized surface of the 
sphere. Since J=(e—1)E/4a 


E. (1 a) 


On the other hand J=N’aF, where N’ is the number of molecules 
per c. c. of the medium, and « — their polarizability. This leads to the 
famous Lorentz formula: 


£c don : 
£4-2 se ue ie (2) 
Since N'—N R (N — Avogadro’s number, M — molecular weight, 
d — density of the medium), 
e—l M 4n 
ee (22) 
The expression 
4 
B= Sy a (2 b) 


is called molar or molecular polarization. 
If the molecules are polar (dipole moment u) then, according to 
Debye’s theory, we have instead of a the expression 


él 
OF 3T 
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and for the molar polarization the expression 


An u? 
OMNE ER, ‘ 
N letz). (3) 
The part of polarization due to dipoles 
; ATN u? 
dip — = ^. | 
3 IKT 2) 


is called „dipole polarization“. 


It is well known that the expression = 1 for a polar liquid 


is much smaller than that for the vapour of the same substance at 
the same temperature. This fact is explained as a result of association 
or coupling of the molecules. Onsager (1936) tries to explain this fact 
in a different way. He supposes the existence of a cavity field different 
from the Lorentz field, and so he obtains for the polarization P a for- 
mula which is different from (2a). 

Onsager imagines the molecule to be a spherical cavity of radius « 
equal to the radius of a molecule in a uniform and continuous medium. 


In the center of this sphere he localises the molecular momemt 5 
(permanent plus induced). Inside this sphere there is a rather weak 
field 

den 


Qu mom 
ND 


(4) 


> 
The moment m, inducing on the surface of the sphere an electric charge, 
forms an additional „reaction field“ 


SS 
> 2(e—1) m 


241 a? G 
so that the total internal field is 

"-—— 

F=G+R. (6) 
If we consider still that 

> as eas 

m= al + u, (7) 
where 

=" — a (8) 

not 2 


(Neo — refractive index of the medium for infinite wave-length), then 
we have cited all Onsager’s assumptions leading to his equations. 
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I do not quote these formulae here, as they will be given in a more 
general form in the next paragraph. Further, Onsager’s formula has 
a rather inconvenient form which makes it impossible to recognize 
the role of the parameters @, R and a in the phenomenon of dielectrie 
polarization. 

The rather serious defect of Lorentz’s theory, namely that it does 
not take into account the action of the nearest molecules, is not re- 
moved by Onsager’s theory, which treats the neighbourhood of each 
given molecule as a continuous.medium, a procedure which, according 
to my opinion, is inadmissible. This is proved by many facts, especially 
by the behaviour of the electric saturation and the Cotton-Mouton 
phenomenon in liquids which will be considered in later parts of this 
work. Therefore, in a generalized Onsager’s theory, the author sur- 
rounds the molecule under consideration with a larger sphere, containing 
several of the molecules in the immediate neighbourghood; the mole- 
cules outside this sphere are treated in Onsager’s way but the few ones 
inside the sphere keep their individuality; they can be coupled in 
their own way causing an alteration in the observed phenomena. 


8 2. Foundation of Onsager’s extended theory. Accor- 
ding to Onsager's theory a polar molecule is placed in the center of 
a hollow sphere of radius a, equal to the radius of the molecule. This 
sphere is contained in the medium which is considered to be contin- 
uous; its surface is to be „smooth“, which seems rather unreasonable 
if we consider that it is built up of molecules as large as.the sphere 
itself. According to Onsager the given molecule induces charges on the 
surface of the sphere, producing a „reaction field“ R. In consequence 
of this „smoothness“ of the surface, the reaction field. does not exert 
any directing action on the molecule, but only increases its permanent 
moment by an additional induced moment. We can say that the re- 
action field merely „stretches“ the dipole, without rotating it. The 
dipole in its rotation would not have any position of minimum potential 
energy. Such a supposition seems to be very improbable and, as will 
be seen later, cannot be reconciled with several facts, especially with 
the behaviour of the Cotton-Mouton molecular constant in liquids. 

Generalizing Onsager's theory, I assume that the sphere sur- 
rounding the molecule considered has a radius larger than a, and thus 
several molecules can now be found inside this sphere. If the radius 
of the sphere grows, the reaction field due to a centrally located molecule 
becomes weaker. But, since the molecules are distributed throughout 
the sphere, the average reaction field per molecule will be stronger 
than that of the central one, and will depend upon the properties of 
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the molecules, and especially upon the position of the moment inside 
the molecule. In equation (5) there appears a factor k, which can be 
slightly larger, equal to or smaller than 1; its value will be chosen to 
satisfy the experimental data. Thus we shall write the equation (5) 
in the form — 


ER > 
m 
=r, (9) 
where 
IAE 
d Apr a (2a) 


On the other hand formula (4) will remain unchanged. In order 
to give to the calculations and to the final formulae a clearer form, 
we shall write the formulae (4) and (8) in the following way 


> > 
G=gE, (10) 
where 
LT (10 a) 
3 2e+ 1 : 
and : 
a= pa’, (11) 
where 
n2,—1 
` ae lla 
p "uu (11 a) 


In the first part of this work we omit purposely the forces that 
couple together the molecule considered with its nearest neighbours 
inside Onsager's cavity, as well as the forces that couple it with the 
discontinuities of its walls. We shall introduce both couplings in the 
second part of this work. In this paper we deduce some equations 
disregarding the coupling; when comparing them with experiments 
we can find some information as to the kind of coupling forces charac- 
teristic of the investigated liquid. 

According to these assumptions the total field acting on a mo- 
lecule, expressed by the formula (6), is 


ES 
[n guy (12) 


Thus, using the formula (7), we obtain for the total moment of the 


molecule the expression 
> 


> > 
m=al*+ u*, (13) 
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> > > -> : jaga 
where F* and u* differ from E and 4 only by numerical coefficients, 
namely 


> > 
ram, (13 a) 
E 
and 
f (13 b) 
u —1—rp u. 
Combining (13) with (12), we obtain a more convenient equation for F: 
P= lm. (14) 
a 
Now we are able to calculate the moment of the forces acting on the dipole? 
> > > i 
MEE ur, (15) 
or 
M =u Esino, (15a) 


> > 
where 0 is the angle between the directions E and yw. Therefore, the 
potential energy of a dipole W = it Mdo will be expressed by the formula 


W= —u*F* cos 0. (16) 
§ 3. Dielectric polarization. (a) Pure liquids. In order to 


calculate the polarization, we must know the component of the moment 
of the molecules in the direction of the external field: 


mg-— aF*-+ u* cos 0. i (17) 
The polarization of 1 c.c. of the medium is then 
J=N'mp=N [aF* + u* cos 6], (18) 


where the double line indicates the statistical mean and N’ — the 
number of molecules in 1c.c. cos 0 is to be calculated in the well 
known way from the Maxwell-Boltzmann theorem: 


5 cos 0 e=WIkT2 sin 0 d0 
1 L 
cos un Re ME (19) 


3 
e—WikTOn sin 0d0 
0 


1 Onsager calculates this moment in a different way; he puts namely 


> > > 
M=Gxm 


and PS nt equations. This method is not correct, because the prt field 


is not à but P, and the dipole, rotating together with the molecule, is not m but JE 
We have called attention to these facts with Prof. K. Zakrzewski in Nature (1939, b). 
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15: : 
where L(x)= coth x Tris Langevin's function of the variable 


ny 
= ET (19 a) 
Taking into account only the first term of the development, i. e. neglect- 
ing the saturation effect, which appears in very strong fields, we 
obtain 


NE iat 
J (^ -— ser (20) 
Since J= (e—1)E/4n we obtain the relation 
ee Ar m H 
—— ——N' —— 9 
| 5 3 (+) 21) 


lf we apply this formula to non-polar liquids (u-0, n?=e) and 

suppose that k—1, we obtain the well-known equation of Lorentz. 

For polar gases and vapours, where e and n? are very nearly 

equal to unity, r and p are very nearly equal to 0 and g#1; the equa- 

tion (21) turns into Debye’s formula: 
e—l 4x 


== —— («+ ai 


D 
3 3 zl: t 


Thus we see that both theories: of Onsager and Debye (eq. (21)) and 
of Lorentz and Debye (eq. (22)) lead to identical results for gases and 
vapours. both theories agree therefore upon the point that the moment 
determined by the method of measuring the dielectric constant of 
gas or vapour is the „real“ moment of the molecule. 

On the other hand, if we apply the Lorentz-Debye equation 
to dipole liquids we obtain generally a smaller value for «. This dimi- 
nution is attributed to association or coupling. Onsager's theory, 
however, suggests that it can also be applied to liquids without any 
further restrictions, and that equation (21), applied to liquids, should 
also give the real value of the dipole moment of the molecule. Prof. 
K. Zakrzewski and the author (1939, a) have already shown that 
this is not so. 


(b) Solutions. We shall consider here solutions of a polar solute 
in a non-polar solvent. The equation (21) applied to such solutions 
changes into 


ml fry lt nl) 085 
3i» = L—Tyo Pı 3kT) 1—?y P2 
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The indices 1, 2 and 12 refer, as usual, to the solvent, the solute and 
the solution. Attention must be called to the fact that the value of 
the moment u* depends not only on the properties of the molecule 
but also on those of the solution, namely on its dielectrie constant, 
since 


u 
à l1 — riz P2 
where 
2(£&12— 1) 
== 24 
Yi NES (24 a) 


N; and N; indicate the corresponding numbers of molecules of both 
components per 1 c. c. of the solution. If c, and e; indicate the weight 
concentrations of both components .M, and .M, their molecular 
weights, dia the density of the solution and N — Avogadro's number, 
then 


Catia « (25) 
Sometimes instead of weight concentrations c, and ¢,, it is 


more convenient to use the „molar fractions“ of both components 
determined by the equations 


Ni N; 


ATMAN? POMENE Ne 
then 
_hM; fo Ms 
ui ie DU e i 
where 
Mi—f M, f M,. (27a) 
Using these equations, we obtain instead of (23) 
ea + 1 Me án h pr? Í 
edere ac EN (a ee 
ES di ER er tz) 1 1— riz P2 Gi 


But Lorentz’s equation applies to nonpolar liquids, and we may 
write it in the form 


— = —Na=P.. (29) 


Besides this we introduce in (28) the expression (3) for molecular pola- 
rization of the solute: 


Pm RN (ut sc py 
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For this purpose the polarization P, must be divided in two parts: 
the non-polar part 
4n 


P=- Na (31) 
and the polar part 
di 4nN u? 
E TM (32) 


The first one is usually determined in the optical way from Lorentz's 
relation by extrapolating n, to infinite wavelengths 


DEDE 
HS 1 M, 


As 2 D 
N90 - 2 d; 


(33) 


It can also be determined by measuring the dielectric constant e, in 
eondition in which the dipoles are immobilized, for instance in soli- 
dified substances 
UBER M, 
" Eg + 2 d, 


(34) 


This method was applied by the author in his previous work (1933). 
The polar part of the polarization, P2?, is determined either in 
gases, according to a formula analogous to (22): 


Ir TH dip Eu M _4n | Le | 


or, more frequently, in solutions. Aecording to Onsager's theory the 


formula (28) should be used for this purpose. 


By using the equations (29), (31) and (32), we can simplify this 
formula and write 


n dip H he 
HPE Ps derer. en (36) 


E — "3 p1 


This form will serve to determine the molecular polarization pj 


of the solute. The quantities p, and p, can be determined from the 
formulae (29) and (33): 


d nd 
puo Ping? P= Page (36 a) 


If we assume the molecular field to be the Lorentz field P = (e+ 2)H#/3, 
then we must insert in the equation (36) the values g,.=(€.+1)/3 and 
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7,,— 0. It is then transformed into the formula derived from the Lo- 
rentz-Debye theory expressing the additivity of polarization: 


Japala Pajas (37) 
where: 
&g—1 Mi, 


£g +2 de 


P= 


(37 a) 


Thus we see that in the Onsager-Debye theory the molar polari- 
zation does not fulfil the additivity relation. 


(e) Interpretation of experimental results. The theories considered 
lead to different ways of caleulating the molecular dipole polarization 
PẸ, but the theory of Onsager supplies different values for pi 
according to the choice of the reaction-field coefficient k. We shall 
make the calculations for mixtures of nitrobenzene and benzene, for 
which the dielectric constant, as well as the other effects of molecular 
orientation, have been carefully investigated. Thus we shall be able 
to apply the same mixture for the calculation of different molar constants 
of nitrobenzene, such as Kerr's constant and others. 

For calculating P2? the author has used the measurement of 
L. Kozłowski (1938). In table I are compiled the values of z and diz, 
interpolated for molar concentrations in round numbers, and pi» 
ealeulated according to the theory of Lorentz (equation (37)) and 
Onsager (equation (36)) for k— 0,5, 1, 1,5. For the non-polar part of 
the molecular polarization P$ the author has accepted the value 41 c. c. 
resulting from his measurements of the dielectric constant of soli- 
dified nitrobenzene (1933). 


Table I. Dipole part of molecular polarization of nitrobenzene as a function o 
concentration (solution nitrobenzene — benzene, temp. 20°C). 


R according to the 


Pob aecording to the theory of 


fe Exe dj N,x10-" eras pase of 
Lorentz t nsager 
i 0, ee, Lorentz ^ 


0 2,282 0,8791 0 


(e. 330) (e.520) (e.240) (c.20) |(e.0,9) (c. 0.68) 
300 472 214 


0,01 2,51 0,883 0,0678 37,7 0,820 0.584 
0,02 2,75 0,887 0,130 294,5 412 212 82,4 0,804 0,579 
0,05 3,40 0,899 0,338 247,4 325 180 82,6 0,676 0,492 
0,10 4,56 0,9174 0,673 207,9 300 159 71,0 0,568 0,434 
0,20 7,05 0,954 1,32 157,5 292 143 58,7 0,430 0,390 
0,30 9,82 0,991 1,95 126,5 293 140 52,6 0,345 0,382 
0,40 12,80 1,025 2,57 105,8 300 138 48,0 0,289 0,376 
0,60 18,80 1,087 3,73 79,0 317 143 47,9 0,216 0,390 
0,80 25,80 1,1467 4,85 63,2 343 151 49,1 0,173 0,412 


1 34,33 1,2032 5,88 


52,8 379 164 51,2 0,144 


0,448 
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These results are illustrated by the curves in fig. 1. Contrary to 
the previsions of Onsager’s theory, the polarization and also the value 
of the polar moment are not constant, and for a pure liquid the former 
is not equal to the polarization in the gaseous state (gas P2"— 366 c. c. 
at t= 20? C). The value of y, determined for the liquid state, will not 
be the true value (see C. J. Bóttcher, 1938). 

We see that both in Onsager's theory and in that of Lorentz 
we are obliged to assume a certain factor which alters the polar orien- 


0 023 OMS! ‚087 =D 


Fig. 1. Dependence of dipole polarization on concentration by the assumption of 
a Lorentz or of an Onsager field. 


tation. According to previous workers this factor was the association, 
i. e. the formation of a certain number of double molecules with a re- 
sultant moment equal to zero. According to Debye (1935) this factor 
-is a quasi-crystalline coupling; we shall call it the type I coupling. 
But many facts have compelled the author (1937 a; 1939) to intro- 
duce a coupling of a second kind (type II coupling), differing 
from the association proper and consisting of a stronger or weaker 
coupling of each molecule with one of its nearest neighbours. Exper- 
imental facts convince us, that the type II coupling is very important, 
though the type I coupling is irreplaceable in some ways (see part II 
of this work). | 
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According to these opinions the interpretation of the results shown 
in fig. 1 would be approxımately as follows: 

Lorentz curve. This curve extrapolated for infinitely great di- 
lution (f,= 0) gives the value cobs? equal to about 330 c. c. It is smaller 
than the value pee and indicates the existence of a coupling of nitro- 
benzene molecules with the solvent. The corresponding values of the 
reducing factor: 

s P2" 
Pai 
can be read on the right hand scale of the diagram. From this value 
the curve begins to fall. The reducing factor 


dip 
er (38) 

Pras 
decreases. It proves increasing coupling of nitrobenzene molecules 
among themselves. We pass over the question of the nature of this 
coupling; in every case it diminishes apparent molecular moment. 
Results to be discussed later, especially the double refraction in a ma- 
gnetie field, show that both kinds of coupling are present. Anticipating 
the results of the second part of this work, we can mention that if we 
adhere to Lorentz’s theory, we shall meet with insurmountable diffi- 
culties. We shall be compelled to reject that theory and to adopt an 

extended Onsager’s theory. 


Onsager’s curves. (1) For k=0,5, oP $” ig equal to about 500 c. c., 
- which proves that the coupling with the solvent raises the average 
molecular moment above its real value; this increase grows as k gets 
. smaller. As such a behaviour seems to be deprived of physical meaning, 
we exclude k equal to 0,5 without further consideration of the curve. 

(2) For k—1 the coupling with the solvent is a little stronger 
than according to the curve of Lorentz. The coupling grows with con- 
centration and the curve drops to a minimum at f equal to 0,4. It 
rises again for higher concentrations, but quite slowly. This seems to 
indicate the existence of an additional coupling of a third kind (type III), 
depending on the formation of aggregates with increased moments. 
It is not supposed that they are quasi-parallel pairs (fig. 2 a), since 
nitrobenzene molecules show inclination to. quasi-antiparallel coupling 
(fig. 2b). We should rather suppose that they form a combination 
of two quasi-antiparallel pairs, coupled quasi-parallelly (fig. 2c). 
Again, it should be mentioned that this opinion is strongly supported 
by the Cotton-Mouton effect in nitrobenzene solutions. We may 
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add that the above mentioned aggregates loose and transitory rather 
than rigid, and as such they will be considered in the theory of coup- 
ling in the second part of this work. 


E cerei e 


a 


Fig. 2. Coupling of two dipoles: 


Te (a) nearly parallel, (b) nearly antiparallel, and 
(c) nearly parallel coupling between two nearly 
antiparallel pairs. 


— 


(3) For k=1,5 the coupling is very strong, as the curve P3? 
runs very low. In particular, as the dilution becomes infinitely great 
(fa— 90), the coupling grows immensely, which again seems to be without 
physical meaning: we exclude therefore the value k— 1,5 in Onsager’s 
theory. : 

In eonsequence of the above discussion, the most suitable value 
for the reaction field coefficient is k— 1 or a value very little different 
from it. In the sequel we shall admit k= 1. In Onsager’s original theory - 
it means that the radius of the spherical cavity is equal to the radius 
of the molecule. But let us remember that in the extended theory 
the sphere must have a larger radius, as it must contain several mole- 
cules. This cireumstance reduces the reaction-field coefficient k, but 
on the other hand another factor tends to increase it as the polar mo- 
ments cannot be considered as small elementary moments acting 
from a great distance; they are not placed exactly in the center of the. 
sphere and may approach closely to the surface of Onsager's cavity, 
and thus produce a much stronger reaction field. 


c 


8 4. Kerr Effect. (a) Pure liquids. Kerr's constant is given by 
the formula TE 


Ne—Nx 1 
E 2 3 
Ku | E: ( 9) 


where n is the refractive index of light and nz and nx are the refractive 
indices for light vibrations in the direction of the vector of external 
electric field E, and in the direction X normal to the latter respectively. 


It prooved to be convenient to introduce in connection with the Kerr 
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effect in particular in solutions, the so called „molar Kerr constant“ Kr. 
Several authors give various formulae for this constant, differing 
by a numerical coefficient. The following definition seems to be most 
suitable 

Eg —HEx 


K™ = yz ; 


(40) 

where Rg and Rx are the molecular refractions in the directions of E 

and X, à SEE 
Kew ecl Me LL e (41) 
FomL-2 d’ n,+2 d" 

and F is the intensity of the local field. 

Taking into account that nz and nx differ only slightly from m, 
we come to the following formula connecting K™ and K 
en? M (=) 


Kem amp m = 


T (42) 
If we assume that the Lorentz formula for the internal field (1a) 
is true, we get the well-known formula 
En 
Ke a Ip & bern: 
On the other hand, in the case of Onságer’s field we define the molar 
Kerr constant K”* by replacing F by the local field which acts mainly 
in orientating the molecules, i. e. by putting @=gE. (cf. the formulae 
(10) and (10a)). Thus, we have 3 


AA Rr —R IR 6n? M1 : 
Se Ara es 

The Kerr constant may be also expressed by means of con- 
stants characterizing the molecule itself (and not the medium, as the 


constants e, n, K). For this purpose we make use of the Lorentz — 
Lorenz function well known in optics: 


(42 a) 


Km 


dispu LT MEC 
Re=7-Ny, and Eı=—-Ny, (41 a) 


where y, and y, are statistical averages of the optical susceptibility 
of the molecule. Thus, for Onsager's field we get 


= NEE. (44) 


It remains now to calculate the susceptibilities and their statistical 
averages. 
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We shall denote by b, bz, bz, the components of the optical suscepti- 
bility of the molecule, referred to its main axes %,y,2. Let the external 
electrie field E act in the direction of the Z-axis of another coordinate 
system X, Y, Z fixed in space (fig. 3). The angles between the molecule 


Fig.3. Explanatory diagram for Kerr effect. X and Z — directions of vibrations 

of the light vector, Y — direction of the beam, Z — direction of the electric field, 

%, Y, z — coordinate system connected with the molecule (z — axis of symmetry 
and of the electric moment), 6, o, y — Euler’s angles. 


axes and the field E are o,,0,,04, and those between the molecule axes 
and the direction X are a4,05,03. We may consider two cases in which 
the electric field of the light wave E, is respectively parallel and per- 
pendicular to E. In both cases an induced electrie moment will arise 
in the molecule. The components of this moment along E and X respec- 
tively are (1) 


for E,|E: 
m,— b, Eoy cos a, - b, Eoy cos a,+ b, E,, COS a, — yg E,, (45) 
where 
yg b, cog a, + b, cos? a, Fb, cos? a. ; (45 a) 
(2) for E,| E: 
m y= b, E, cos a, + b, Eoy cos a, + b, Eo C08 a; — y 4 E,, (46) 
where 
y y= b, cos? a+b, cos? a+, cos? ay. (46 a) 


Let 0, p, v, be Euler’s angles determining the position of the mole- 
cule in reference to the system X YZ (as shown in fig. 3); 0 and 9 cor- 
respond to geographical longitude and latitude, y is the angle of rotation 
around the z axis. We have 

cos a,— — sin 0 cosy, cos a—sin Ü sin y, COS a5— COS 0, 
cos aj= COS 0 cos p cos y—sin o Sin y, (47) - 


C08 a$— —cos 0 cosy sin y—sin p COBY, COS as—sin 0 cos p. 
| m 
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We shall assume that the molecules are rotationally symmetrical 
and that z is their axis of symmetry. Thus, the rotation of the molecule 
around the z axis is completely free in spite of the external field E, 
and all possible values of y are equally probable. We may therefore 
by inserting the direction cosines (47) into (45a) and (46a) calculate 
the geometric mean over all positions around the axis 2: 

270 1 27 
z; | vet and z; frw 
0 0 
However, all positions of the z axis relative to the direction of the 
electric field E also have the same probability. Thus, we calculate 
the mean value for all values of the azimuth o by means of similar 
integrals, but in reference to gp. Putting b,— b, we get 


yg b tH (b, —b,) cos? 0, y,-—4$(b,--b,)— $(b,— b) cos? 0. — (48) 


Finally, the statistical means of y, and y, should be caleulated by 
calculating the statistical means of cos? 0 over all possible positions 
of the axis of symmetry z, which now are obviously not equally probable, 


D b,+(b,—b,) cos? 6, y 1 (b, 4- b,)—4(b,—b,) cos? 0. (49) 
Henee, thanks to (44), we have 


Km* — InN 3 


bs Gi = 1 (cos? 0 0— 1). (50) 


Thus, we have reduced the problem to the calculation of cos? 0. The 
Maxwell and Boltzmann theorem yields 


f cos? 0 - e-Wir sin 046 
cos? I. (51) 
fe sin 646 
0 
To calculate this expression the potential energy W of the mole 
must be calculated. l 
Let the components of the polarizability of the molecule in an 
external field be 4,,0,,0,; the molecule being symmetrical around 
the z axis, we have a,—4,, u,— 45-0 and u4— p. If the field F 


acts on the molecule, the moment of the morgua will have the com- 
ponents 


m.—a, Fx, my=q4F,, m,—a,F, t p. (52) 
According to Onsager's theory, the molecule should be treated as if 
enclosed in an empty ellipsoidal space of identical dimensions with 
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those of the molecule. In the modified form of this theory, however, 
we assume that several molecules are found in a spherical volume the 
radius of which greatly exceed that of the molecule. The internal field 
is expressed as before by (12): 

3 r r % 

FE. m. Fy=gE,+ 3my, F,—gE.Yom. (53) 
The solution of the equations (52) and (53) may be presented in vector 


> E 
form, where the vectors F and m will be composed of two components: 


> > 
one parallel to EZ, and another parallel to u. We shall use the following 
symbols 


Pe—.—9 [n (54) 
1—r79A,.! 
x: T > 
le ue (55) 
where 
a=, =% (a-iimd ara). (56) 


A, and A, give a measure of molecular anisotropy. 
Further let us put 


ag = xaz 
with À 
l—T"pA 
LE N T 


> > m m 
Finally, let u stand for a unit vector parallel to u. The solution is as 
follows 


mem leat (ya 60s Js (58) 
3 

goriiss > 

m — a, F* 4- [u*+ (af —a,) F* cos O]u. ! (59) 


For isotropic molecules (a= 4,— 43— a) the above formulae become 
identical with our previous formulae (14) and (13). 
The moment acting on the molecule is 


> > > i 

M=Fx [u*+ (aš —aı)F* cos 0]v, ~ (60) 

therefore l 
M = u*F* sin 0+ (a$ —a,)F* cos 0 sin 0. (60 a) 

Hence the energy 


* 
. W=- f Mao= wir" cos 0 - Zr cos? 0. (61) 
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For isotropie molecules this expression goes over into (16). The form 
of (61) is such that without the asterisks we obtain the same expression 
as with the Lorentz field; in this case F means of course (e+ 2) E[3. 
We may now calculate cos? 0 from (51). Expanding the exponential 
into a series and integrating term by term, we obtain 


—— gl 2 4 

cog? 0— 3 T 4g 7" o 3g 9T - (62) 
where ER A 

at, = pe, (63) 


From the above and from (50) the final expression for Kerr’s molar 
constant is obtained. We give it here in symbols which are connected - 
with the theory of the Kerr effect based on the Lorentz field: 


1 
m* — © * * 
K nA PLE Ple rta (64) 
where 

2 
f= tap $) (b —by), | 

i (65) 
02— ame” (53 bi) | 


In the case of a Lorentz field, i. e. if we put g= (e+ 2)/3 and r=0, 
the asterisks in the above expressions should be omitted and we obtain 
the well known formula 


K™ = 2xN(6, - 6,), (66) 
in which the expressions 
x 2 
0,— 3g p 6892) (b,—b,) 
and 
0 2 2 
2— gp pga (0s —b1) (67) 


are a particular case of the expressions known from the theory of the 
Kerr effect, namely 


1 
(,— ggg T a) (b, —b2)+ (a — a5) (ba —bz) + (a3 — a) (b —b,)] 
and | (67 a) 
1 à ~ 
b= +. 
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(b) Solutions. Kerr's molar constant Km* is calculated in the 
previous paragraphs on the one hand, by means of quantities characte- 
ristic of the medium (formula (43)) and on the other, by means of 
quantities characteristic of the molecule (formula (64)). In the latter 
formula, however, there appears the quantity r which is characteristie 
of the medium, that is the reason why the constant K™* is not additive 
for solutions, in opposition to Km which is based on the assumption 
of a Lorentz field. Nevertheless, we shall calculate K”* for solutions. 
When taking the solution to be a continuous medium, we shall use 
formula (43) but, regarding it as a mixture of two kinds of molecules, 
we shall use formula (64). 

We shall limit our reasoning to the case of a polar substance 
in a non-polar solvent, the anisotropic term 6, of the polar substance 
being so small in relation to the term 0, as to be negligible. A similar 
reasoning to that of $3 (b) leads to the application of the formulae 
for K”* to the case of one gram molecule of the mixture (M, for- 
mula (27a)). We obtain 


nsu qM. 2n N0T 2z N06, 
ae a a Ea ee 
(niat 2) giz diy — (1—T,339,4,) (1— ryp3 4, * (1— 7 p, 4) 
As the first term on the right-hand side is a small fraction of the second 
for both big and small concentrations, we may put in its place the 
expression 


2z.N 0, j 
C 
Regarding 
2znN0,— Ki | (69) 
and 
2nN0,— Ko’, (69 b) 


as Kerr's molar constants of the dissolved substance and the solvent 
respectively, for a Lorentz field, the formula (68) may be written in 
a simplified form as follows 


K” Je) 


LH MEE Bleed MEUS EEE gah c uos 10 
En (L—713p1)? ht (rap? np. AD fz (10) 


where j 
6n?, 1 M fs 


* (n3, + 2)" 9%, d, ’ 


ee dr MS 


Ki =K, (10 a) 


56 A. Piekara 


By means of this formula we may calculate K5', or 2zN6,. Accor- 
ding to Onsager’s idea this quantity (as well as the dipole moment) 
should be invariant for all concentrations and for vapour. We shall 
see that this is not so. 

We notice that in the case of a Lorentz field, i. e. if we put ri = 0 
and g,,—(454-2)/3, formula (70) assumes the well-known form 


Kn=Kıfhıt Krk, (71) 
by which the additivity of Kerr’s molar constant is expressed. 


(e) Interpretation of experimental results. We shall apply the above 
results to the case of a mixture of nitrobenzene and benzene. Nitro- 
benzene actually fulfils the condition assumed in (b), namely? 6f < 07. 
Thus the formula (64) takes the form 


205N 0,74 Wee Km 
(1—rpA) (1—rpA)* (1—rpd,)? (1—rpÀA)* 


Hence. we may calculate Kerr's molar constant K” for pure nitro- 
benzene. However, Kerr's molar constant of nitrobenzene determined 
from solutions is to be calculated (as Kz) from (70). In the case where 
fe=1 this formula obviously takes the form (72). 

We shall take numerical data from experiments of H. Friedrich 
(1937) performed at 23? C. Friedrich measured Kerr's constant RK, 
of nitrobenzene-benzene solutions and caleulated the molar constant 
of the solution „K,“ from a formula which differs from our formula (42a) 
only by a factor of 9, namely $,—41 Kiz. Accordingly, the constant Kia” 
appearing in (70) is calculated as follows: A5*—$,(e,,-- 2)?/g?, ; its 
value is tabulated in the sixth column of Table II. Kerr's molar con- 
stant for pure benzene K7 is 98, for f, — 0 and amounts to 48,33 x 10-12; 
p, and p, are calculated from (36a). For P? we have assumed the polar- 
isation of solid nitrobenzene, i. e. 41 c. c. Finally; we have evaluated 
Ay and Ag on the basis of b,, b, and b, taken from Stuart’s book quoted 
previously?. Table II contains the results of calculations of Kerr’s 


fmt — (72) 


* According to data from investigations of the Kerr effect in nitrobenzene 
vapour at 235,5? C, we have 6,/0,— K,/K,—7/139— 0,050 (see H. Stuart Molekül- 
struktur, Berlin 1934, table 55, p. 211). At 23° C this ratio is 0,029. The ratio o*/ox is 
still smaller. 

— 3 In table 59, p. 222 we find b,—132,5x 10-5, b „=77,5x 107” and b,=177,6x 1075. 
These are polarizabilities of the nitrobenzene molecule for visible light determined by 
means of the Kerr effect and the depolarization of light in nitrobenzene vapour. As- . 
suming an approximate proportionality between a,,qa,,a, and bi, b2, b3, we calculate 


b 
TA 4 (ba + ba) —0,813 and A7/— by 


TUE b, + bs) 3 1(b,+ b. ee chos 
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molar constant of nitrobenzene according to both the generalised 
theory of Onsager (for k=1) and to the theory of Lorentz. In the 
first place we see that the „constant“ Kz’ is not constant as it should 
be according to Onsager’s theory but depends in a rather queer way 
on the concentration. This dependence is represented by the diagram, 
fig. 4. 


Sa 
Hl! 
| DR 
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Fig. 4. Kerr's molar constant of nitrobenzene as a function of concentration. 


Table II. Kerr's molar constant K7' of nitrobenzene for different concentrations (from 
nitrobenzene — benzene solutions at 23°C). 


Onsager Lorentz 


; , = * 9 
Ip E dà NXI0 Ryyx 1008 Ky x10 K™10 R [Kao R 


0 2,26 0,8748 
0,0930 4,42 0,9097 
0,1957 7,01 0,9468 
0,2618 8,80 0,9701 
0,3532 11,3 1,0014 
0,4659 14,6 1,0387 
0,5769 18,3 1,0743 
0,6872 22,5 1,1083 
0,7829 25,6 1,1369 
0,8805 29,6 1,1654 
1. 34,5 1,1992 


* Friedrich’s data multiplied by 9. 
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The reducing factor R is a measure of the diminution of the 
molar constant in comparison with the value for vapour. Ris calculated as 
the ratio of the constant K7 determined from the solution or from 
the pure liquid to the molar constant Toros of the same substance 
in gas or vapour form: 

ES 


MEI 
K gas 


(73) 


For nitrobenzene vapour at 239 € K5,—11,736 x 10. The calculated 
values of E may be found in Table II; the same diagrams on fig. 4 
represent the dependence of R on f,. 

The Lorentz curve, the form of which is well known, shows that 
„the molar constant of nitrobenzene falls abruptly with increasing 
concentration, so that for pure liquid the reduction factor is only 0,054. 
It might be noticed that the decrease of the value of K3' is not connected 
with the decrease of any quantity directly observable by the experi- 
menter. Rather the opposite occurs; the experimenter when measuring 
the ordinary Kerr constant states that it increases remarkably with 
increasing concentration. Only when he calculates the molecular con- 
stant fu, by means of a formula analogous to (42a) does he divide 
by the term (e+ 2)? which increases strongly with increasing concen- 
tration. This effect is still greater in the calculation of electric satura- 
tion ($ 5), where the divisor is (e+ 2). | 

It might be said that this remarkable fall of K5' is due to the use of 
the Lorentz field, which is expressed in the formula (42a) in the form 
of a very large divisor. One is inclined to arrive at the conclusion that, 
instead of explaining the fall of K3" by a very strong coupling of mo- 
lecules 5, it would be more natural to revise the method of calculating 
the internal field. In fact, assuming the field of Onsager, one obtains 
a fall of the curve K?' which is much less steep, and though some degree 
of coupling must be aecepted as an explanation the coupling may be 
much weaker. 

Nevertheless the Onsager curve has a very remarkable form; 
it corresponds exactly to the variation of the molecular polarisation 
(fig. 1). At first the curve falls down, which is an argument for an 
increasing antiparallel coupling of the dipole molecules of nitrobenzene. 
For stronger concentrations there is an overwhelming influence of 
another type of coupling which, I am inclined to suppose, consists 
in a parallel coupling of antiparallel pairs (fig. 2). 


* According to Friedrich, this value is one ninth as large, i. e. 1,301x 107°. 
5 As has been done by Debye (1935), Friedrich (1937) and the author (1939). 
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An interesting irregularity appears in the neighbourhood of 
f=0,3—0,4, which is connected with the attaining of a maximum 
by Kız, found experimentally. 


$ 5. Electrie Saturation. (a) Pure Liquids. The component 
of the moment of the molecule along E is given by the formula derived 
directly from (59): 


Mp — [a,4- (at —a,) cos? 0] F*+ u* cos 0. (74) 


Our problem consists in caleulating the statistic mean value of 
mpg retaining the terms with F°*. For this purpose we shall calculate 
cos Ü and cos? 6. The latter value has been calculated already and 
is given by (62) and (63). We shall calculate now cos 0 using formula (61) 
for the energy of the molecule. With the same meaning of the symbols 

n (63), we have 


FLA 
* — YF 2 ^ 
f «os DE 0 212 ne 
1 4 
cos a ae et 
er k ES * 2 d 
fe x" COS 0 — X%7 COS ? sin odo 


1 * 
ago f.. (08) 


By substituting the expressions for cos 6 and cos*6 into (74) 
we get 


ms (a4) "+14 20809), (76) 
where 
24 + az gt .. 2(a4$ —aJ* ] 
tvm sanae Urge CIT UNES 
TET a is y icut M. 
ES pen oem 
The dielectric constant is expressed by the formula 
eiaa N IE (18) 


Differentiating (76) with regard to E and substituting into (78), 
we obtain an expression for the dielectric constant in a strong 


+124 N' (01 +20% — 03) — ——35 


(79) 
(1 See! a 
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whereas in a very weak electric field we have 


u? 1 | Jo 


3kT (1—rg pA? | 1 —r9 ph 22 


co= 14- 4x N* & 


(the index , refering to the dielectric constant in a weak field). Under 
the influence of a strong applied field the dielectric constant is changed 
by de=e—e,. Using the definition 


1 g 
7’ Re Er pene” 81 
en let an eee (81) 
we obtain by subtracting (80) from (79) 


g? 
As=f(é,...)—f(e,..-) + 12N" (65+ 202 — 08) irme 


For extremely small Ae (and only such are to be taken into account) 
we may write 


2 
fle) fe.) 2: Ae, 
hence 
Á MÀ '(* o0*  px* rie gi d 2 
A= A 12zN'(0f-4- 202 — 03) reves 3 (82) 
de 


In the case of a Lorentz field (g=(e+2)/3, r— 0), the above formula 
takes the well-known form 


4 
Ae 12zN' (0,4- 20,--0,) (5) E, (83) 


where 0, 0, and 0, stand for expressions identical with (77) but with- 
out asterisks. 

Formula (82) assumes a much simpler form for substances with 
a relatively big moment u, for which a is small in comparison to the 
dipole part of the polarizability. For these substances the terms 6, 
and 0, may be neglected in comparison with 6, (in a still greater degree 
this applies to the expression with the asterisk). E. g., for nitrobenzene 
at 20°C we have 6,=0,013x10- and 0,—0,44x 10—?, whereas 
0,— 10,76 x 10-2. 

In these conditions 


dicm 12zAN'0, g? 


E pu Orpa a rph (84) 


where 9f/de is the derivative of (81). 
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If according to our assumption a is neglected, the formula 
for 9f/[9e becomes 


ej ei fl A, 24s 6p 
de Bet] It rr ls (Bes) 


(b) Solutions. We shall apply (84) together with (84a) to the 
solutions of nitrobenzene in benzene. The error introduced by the 
previously mentioned simplification does not exceed 10°/, for diluted 
solutions and decreases almost to 1°/, for concentrated solutions. 
More exact calculations are not needed as measurements of electric 
saturation are unusually difficult and their precision is limited. In’ 
the literature we find only few investigations concerning electrie sat- 
uration. Besides Herweg's fundamental work we have the paper of 
Kautsch (1928) who measured Ae with great care and precision; 
his investigations, however, concern ether and chloroform only. More 
dipole liquids were investigated by Malsh (1928, 8) and Gundermann 
(1930) but these authorsdid not take all the necessary precautions and their 
method gives only qualitative and relative results. Measurements 
of Ae for solutions of nitrobenzene in benzene and also in other solvents 
were made by the author together with B. Piekara (1936)9. The 
results of these measurements will be discussed here in connection 
with theory. 

Formula (84) may be used here in a remarkably simple way, 
Ae for the pure solvent (benzene) being so small as to remain within 
the limits of experimental error; therefore we may write 4e= 0. Con- 
sequently, it is sufficient to replace N’ by N; from (25) to obtain the 
formula for mixtures: 


Sm 
(85) 


ð 
=) (1— ri pA? (1— ri pÀ4)* 


where M, p, 4, and A, refer to the nitrobenzene molecule and o//de is 
expressed by the same formula (84a) with the only difference that c, 
and 7,, should take the place of e and r respectively. The new symbol S7 
appearing here is defined by the formula 


Sm— — AnN6, | (86) 


6 These investigations were continued in the following years and partly published 
(1937). The remaining notes were lost during the war. The measurements with 
pure nitrobenzene were repeated in 1939 by the author and A. Lempicki by means 
of a method different from that hitherto applied and using a cathode-ray oscillograph. 
The outbreak of the war interrupted this work and the apparatus has been destroyed. 
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and plays the same part as K"—=2n.N6, in the theory of the Kerr effect. 
We shall refer to it as „the molar constant of electrie saturation“. 

It is easy to prove that if a Lorentz field is assumed and electro- 
striction is neglected 8” is the increase of molecular polarization per 
unit of external field F: 


Pr P 
Se oF (87) 
Pg is the polarization in the field, P, — without a field: 
— —1 M 
Tus £—1 M Eo (88) 


EUA re 


Assuming Ae= e—e, to be small in comparison with e and applying (83), 
we get 
Sm— Az N (0,4- 20; — 04), (89) 


which for a particular case yields (86). 

We shall calculate $7 for nitrobenzene. Putting u= 4,23 x 10-18, 
we obtain 0,— u*/40k? T3 — 10,76 x 1022, hence S7 — —81,4 x 10—10. This 
value is obtained on the assumption of an absolute mutual independence 
of the molecules. It may be therefore applied only to gas or vapour; 
we shall use the symbol Sas for it: 


D — 81,4 x 10 15, 


If we insert this value into (85) and calculate As, we get a result 
largely differing from the experimental values of Ae... This shows 
directly that Onsager's theory cannot even approximately deal with 
the electric saturation in liquids (no more than with the dielectrie 
polarization and the Kerr effect). This, however, is no reason to condemn 
Onsager's theory as false. We shall consider it as insufficient and we 
Shall try to explain the discrepancies between theory and expe- 
rimental results by the same coupling of molecules as in the case of the 
polarization and the Kerr effect. 

We shall not calculate the values of 4e, based on Sz, but starting 
from experimental values of 42, we shall calculate the value of S= by 
means of (85), and hence the reduction factor R, of saturation: 


ae (90) 


R= 


The more the latter differs from unity, the greater is the part played 
by the coupling. 
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(e) Interpretation of experimental results. Table III contains the 
following experimental data: the weight concentration c, of nitroben- 
zene, the number of nitrobenzene molecules N; per c. c. of the solution, 
the density and the dielectric constant of the solution, and the variation 
of the dielectric constant Ac under the influence of an electric field 
of 70 kV/cm. The last four columns contain the molar constant Sm of 
electric saturation of nitrobenzene and the reduction factor R, calcu- 
lated on the assumption of a Lorentz field 


9 A196, m (£ re 


as well as of an Onsager field (formula (85)). 

The dependence of S" and R, on Na based on the assumption 
of a Lorentz or an Onsager field is represented by diagrams in fig. 5. 
The Lorentz curve shows that the decrease of the molar constant 
of saturation with increasing concentration is here still more strongly 


Fig. 5. Molar constant of electric saturation of nitrobenzene as a function of con- 
centration. A part of the Lorentz curve is presented on a larger scale; the positive 
effect of saturation may be seen. 
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Table III. Molar constant of eleotric $aturation as a function of concentration (from 
solutions of nitrobenzene — benzene at 21° C and E=70kV [em) 


Onsager 


Sm y 1019 E 


Lorentz 


n —21 104 
cĉ — N$x10 diy Erg AEX gm y ire R, 


0,89: 3,000 — T, — 41,6 +0,51 —7,48 +0,092 
0,0984 0,435 0,904: 3,75 — 2,71 —17,0 4-0,21 ' —3,65  --0,045 
0,1307 0,582 0,913 4,21 — 4,06 —14,1 +0,17 —3,36 -+0,041 
0,2874 1,343 0,955 7,20 .— 9,47 — 2,93 + 0,036 —2,11 +0,026 


0,522 2,63 1,028 13,05  —20,3 — 0,45 + 0,0055 —1,67 +0,0205 


0,802 4,41 1,125 23,0 —32,7 —0,057 +0,00070 | —1,35 +0,0165 
0,902 5,14 1,164 28,3 — 15,7 —0,0108 +0,00013 | —0,53 +0,0065 
0,972 5,57 1,192 32,5 +32,7 +0,0124 —0,00015 +0,99 —0,012 
1 5,88 1,203 34,33 +75,4 +0,0220 —0,00027 +2,12 —0,026 


marked than in the case of the Kerr effect (fig. 4). The decrease is 
so abrupt that the section corresponding to higher concentrations 
had to be drawn on a bigger scale. As has been previously mentioned 
in § 4 (c), this is caused by the fact that the ratio F/E=(e-+2)/3 enters 
‘to the fourth power. If we assume an Onsager field, there appears 
another, much smaller, factor (see (85)); hence the decrease of S™ 
and R, with increasing concentration is here much less accentuated 
and in general the reduction factor is greater as shown by the curve 
of Onsager. Thus the coupling of the nitrobenzene molecules is not 
so strong a8 was supposed on the basis of a Lorentz field. This agrees 
with our conclusion drawn from the plots of Kerr's molar constant 
and of the molecular polarization. However, the final course of the 
curve of the electric saturation is different, S changing its sign from 
negative to positive. In the second part of this work we shall discuss 
this effect which is very characteristic of the second-kind coupling 


$6 The Cotton — Mouton Effect. (a) Pure liquids. The 
Cotton-Mouton effect is the magnetic counterpart of the electric 
Kerr effect. Plane polarized light passes through the liquid placed — 
in a magnetic field, the lines of force running in a direction perpendic- 
ular to the light beam, and at 45° to the plane of the vibrations of the 
light wave. The effect consists in the change of the plane polarized 
beam into an elliptically polarized one. The Cotton-Mouton constant 
As eee by the formula 
Wg—nx 1 
e ue iyi (91) 
in"which n means the ordinary refractive WES for the given wave- 
lenght, ny and ny are the refractive indices for. vibrations parallel | 
and perpendicular to the magnetic field respectively. 
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We shall define Cotton-Mouton molar constant, according to 

Kónig (1938), as the difference between molecular refraction for 
parallel vibrations and that for vibrations perpendicular to a magnetic 
field of 1 oersted: 
Ra Rx 
a 
where Ry and Rx are expressed by formulae similar to (41). In this 
case we do not make any distinction between the external and the 
internal field, as the magnetic permeability of most liquids is practically 
the same as that of the vacuum. The formula for C" may be deduced 
from (42) directly, yielding 


om (92) 


2 
ee (93) 


(+22 d 


We see that in this case we are not concerned with the Lorentz 
or Onsager method of calculating the local field. Nevertheless, we shall 
consider the magnetic birefringence, as it will supply us with important 
information concerning molecular coupling in polar liquids (cf. part II 
of this work). The exceptional character of the Cotton-Mouton 

effect in contrast to the other effects discussed is connected with the 
fact that in this case the dipoles play no other part than that involved 
in the mutual coupling of molecules. There is no dipole orientation 
in an external field. Further, from the point of view of experimental 
methods, this effect is much more convenient than the preceding one, 
- especially because electric conductivity does not intervene. 

We shall express the constant C7 yet in another way, namely 
by means of parameters characterizing the molecule itself. We have 
a ready formula for it in § 4, formula (66), where in the case of magnetic 

. dipoles we have to put 0,—0, as the dipole liquids are diamagnetic. 


Thus, we get 


2 
= ipi Te) (b —b). (95) 


"where 


c, and e, are here the magnetic polarizabilities for directions parallel. 
and perpendiculat. to the axis od symmetry of the molecule respec- 
tively. j 

Formula (94) is based on the assumption of complete freedom 
of the molecules and is therefore applicable only.to gases and those 
liquids in which there is no limitation of freedom of the molecules. 
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As we shall see, remarkable discrepancies will appear in cases in which 
coupling. takes place. 

(b) Solutions. As before, we shall limit ourselves to solutions 
of polar liquids (index 2) in a non-polar solvent (index 1). 

From $4(b) we shall take the formula 


Ci — €t fı t O2 fas (96) 
where . a 
6n 
Cr = C 773 = d i 
(niy--2) dy 
and 
m 6n? M, 
Q1 0 


Cz should be determined from the equation (96) and compared with 
the value of the same constant for gas (Cas). Therefore, we shall cal- 
culate the ,reducing factor" Roy: 

Rew= ue (97) 


gas 


(c) The interpretation of experimental results. As far as we know, 
König (1938) was the first to calculate the Cotton-Mouton molar 
constant. His calculations were based on measurements with solutions 
of nitrobenzene in carbon tetrachloride performed by A. Goldet and 
the author with the large electromagnet at Bellevue in 1934. The 
results are presented in Table IV. The first column contains the weight 
concentrations, the second — the densities of the solutions extrapolated 
from measurements of N. N. Pal*, the third — the number of nitro- 
benzene molecules N3 per c.c., calculated from the preceding, the 
fourth — the Cotton-Mouton molar constant as given by König ® 
and the fifth — the reducing factor Rey. In order to calculate the 
latter, the value of the. Cotton-Mouton molar constant for vapour 
is needed. Kónig made an attempt to measure it, but found that the 
magnitude of the effect is within the range of experimental errors. 
For 20? C and 670 mm Hg he obtained 


Ops < 550X107", 


whereas for an infinitely diluted solution of nitrobenzene we obtain 
by-extrapolation to e= 0 


QU = ca 400 x 10^" 


7 From Landolt,Bornstein Tables, vol. III, part 3, p. 1947 (1936). 
8 There is only a diagram in König’s paper. The numerical values of the fourth —— 
column were obtained by careful measurement of the co-ordinates. | | 
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We shall use this value instead of Czas for calculating Rey, the obvious 
result being that Rcy=1 for c,— 0, which is equivalent to the assump- 
tion that the coupling of the nitrobenzene molecules with the solvent 
does not influence the constant C$?'. In fact, as will be shown in Part 11 
of this work, this influence is negligible. 


Table IV. Cotton- Mouton molar constant Cy as a function of concentration (from 
nitrobenzene — CCl, solutions at 209 C) 


d N; x 10721 Ons Mo —— 78 


2 12 CM 
0 1,5940 0 (ca. 400) 1 
0,1217 1,5365 0,914 426 1,065 
0,3084 1,4530 2.191 490 1,225 
0,5461 1,3550 3,619 572 1,43 
1 1,2035 5,885 700 1,75 


As may be seen from (94) the Cotton-Mouton molar constant 
should really be independent of the concentration of the molecules, 
whereas, as shown in Table IV and graphically in fig. 6, Cz’ increases 


EEE Er 


Fig. 6. Cotton-Mouton molar constant as a function of concentration. 
J 


rapidly with the concentration. The „reducing factor“ Rom is in this 
case greater than unity, in contrast to the preceding effects. This 
unusual behaviour is caused by the type II coupling, and the effect 

5* 
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is readily understood from the following exaggerated example. Let 
us suppose. that the nitrobenezene molecules are so strongly coupled 
two by two that one may consider one gram molecule of the liquid 
as composed of N/2 associated pairs. The magnetie and optical polari- 
zabilities of each pair are doubled (if we neglect mutual inductance 
influence). In consequence the expression 0, (95) is increased by a factor 
of four and therefore the molar constant C™ is doubled, hence Ray = 2. 
In reality for nitrobenzene Roy=1,75, which means that the mole- 
cules are not associated in rigid pairs. The theory of the coupling (see 
Part II) gives for Rey a value slightly smaller than 1,75 if we admit 
a coupling energy required by the three preceding effects on the assump- 
tion of a Lorentz field. To obtain the value Rcu=1,75 we must either 
admit a greater coupling energy, which would disagree with the trend 
of the preceding effects, or admit the possibility of à coupling of mole- 
cules in larger aggregates than in pairs. The latter possibility is in 
accord with the variation of the dielectric polarization and the Kerr 
effect, but only if the Onsager field is assumed. We see that although 
the local field does not appear directly in the Cotton-Mouton effect, 
the latter may supply us with valuable information on this subject. 
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ON SOME PHOTON RADIATION OBSERVED BY MEA- 
SUREMENTS OF COSMIC RAYS AT GREAT DEPTHS 
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The problem of great depths 


As concerns the problem of cosmic radiation reaching great 
depths (some hundred metres water equivalent) it is assumed by most 
physieists that this radiation is a meson radiation (perhaps with some 
admixture of protons) accompanied by knock-on electrons and cas- 
cade showers started by these (Wilson 1943). Quite different views on 
this subject are given by Barnóthy and Forró (1948) who assume 
that at depths greater than 500 m w.-e.itis mainly non-ionizing radia- 
tion produced by the decay of mesons which penetrates there. They 
consider that our instruments at those depths are.operated by second- 
aries of that non-ionizing radiation. 

These differences in views are derived from some anomalies 
‚observed by measurements of cosmie radiation at great depths. The 
best known anomaly is that giving the relation between the intensity 
of cosmic radiation and the depth. If we take the integral spectrum . 
of mesons as given by the relation fen 


I=cE’ 


and we also take the losses of energy of mesons to be proportional 
to the depth H which they have traversed we obtain the relation 


I=ecH", 


in which y is equal to about 1,8. In a diagram log J vs. log # we should 
get a straight line with a slope y, whereas if we take into account the 
experimental results, we get a line broken at about 300 m w. e. (Clay 
1939). Up to 300 m w, e. we have aline with a slope of 1,8 and beyond 
this depth the slope changes from the value 1,8 to 2,8. Because of the 
fact that only mesons with energies higher than 10! eV can penetrate 
to this depth, we might conclude that besides the normal loss of energy 


=> 
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through ionization the particles concerned suffer an additional loss 
through some other mechanism. The other explanation might be that 
the spectrum of energy of mesons above energies of 10% eV drops 
quicker than may be expected from the form ber. 

Lyons (1941) considers an additional loss of energy ot very 
energetic mesons by radiation as a cause of this change of the dependence 
of intensity on depth. The mesons of high energies might radiate most 
of their energy in one collision and therefore we should get in the depend- 
ence of the intensity of cosmic rays vs. depth an additional exponential 
factor giving a quicker drop of the intensity with the depth. A closer 
analysis of this process proves that the losses of energy by radiation 


suffered by mesons of energies sufficient to penetrate such depths 


(104 eV) are too small to explain this phenomenon. 

There are now in the literature of the subject two opinions on 
this problem, based on similar assumptions of meson decay, given 
by J. Barnóthy and M. Forró (1948) and K. Greisen (1948). Bar- 
nothy had already considered this question in 1940. 

The mean range of a meson in vacuum (because of the decay) 
is given by the formula 
K= T° = 
where r, is the mean proper lifetime of the meson, P its momentum 
and uits rest mass. In the spectrum of mesons we have a lack of particles 
of low energies, but in the spectrum of decay products we have a lack 
of particles of very high energies. This last phenomenon will be observed 
for those energies of mesons for which the mean range is equal to the 
height of production of mesons in the atmosphere. In the spectrum 
of decay products of mesons we have to add to the form E !? the 
factor E^! for energies corresponding to the decay range equal to the 
height of the meson production layer. If we assume that at great depths 
we register the products of meson decay and that the losses of their 
energies are proportional to the depth they traverse, we get for suffi- 
ciently high energies a dependence of the intensity on the depth in the 
form E ?* 

The main difference between Barnóthy and Forró and Greisen 
is that Barnóthy and Forró assume that at great depths we observe 
neutral products of decay of normal mesons (ma 200); Greisen on 
the contrary introduces the decay of mesons of mass 300 with a much 
lower lifetime (6,107? sec) and obtains in this way for the spectrum 
of energy a change of the form E !? to E-?? for energies greater than 
10" eV. This is the energy necessary for a meson to traverse some 
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hundred metres water equivalent. Thus, at great depths we observe 
normal mesons which themselves are already decay products. 

For the evaluation of this critical energy, ER assumed 
to long a lifetime for the meson of mass 300, viz. 6.107? sec instead 
of about 8.107? sec, as it is now known. For thè latter value for v, the 
change of the exponent in the energy spectrum should occur at about 
10% eV rather than at 104 eV. We might thus take the explanation 
of Greisen as quite satisfactory, because of the great inaccuracies in 
the data we have here at our disposal, arid because of the neglect of 
the losses of energy due to radiation !). 

Thus, we can assume that in cosmic rays at great depths we 
have to do with secondary particles produced by the same processes 
which are known from experiments at sea level. 

But Barnóthy and Forró in the paper mentioned (1948) 
report among many interesting results also a phenomenon which they 
inerpret as an argument for the opinion that at great depths we have 
to do especially with neutral meson-decay products. This phenomenon 
mentioned already in 1939, was also communieated at the Cracow 
Cosmie-Ray Conference in Oct. 1947. They observed at the depth 
of 1000 m w.e. that the radiation there measured was able to give 
many times more double than threefold coincidences in the same 
telescope. The ratio of threefold to double coincidences in their measure- 
ments was 0,05. Barnóthy and Forró express the opinion that this 
phenomenon proves the existence at great depths of a scarcely ionizing 
component of cosmic radiation. For this component the probability 
of threefold coincidences ought to be decidedly smaller than that of 
double ones. 

The aim of this work, of which we are giving proyisional results 
in this paper, was to check this interesting phenomenon and inve- 
stigate more closely the properties of the radiation which produces it. 


Description of the aparatus 


_ Our measurements were performed in the salt mine at Wieliezka 
near Cracow. The deepest level at our disposal was 282 m below ground. 
The layers above the telescope consisted of rock-salt, sandstone, gray 


3 1 Sinee this paper was submitted for publication, we came across a paper 
of Satio Hayakawa and Sin-Itiro Tomonaga, Phys. Rev., 75, 1958 (1949), 
in which these authors, considering the decay of mesons of mass from 287 to 217 
with a lifetime of 1:10—9 sec and taking into account the losses of energy due to © 
. radiation, have obtained a curve of the dependence of wa. cee! vs. depth in accord- 
f ance with experimental data. 
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silt and silt with grains of rock-salt, gypsum and anhydrite. The water 
equivalent was estimated from samples taken from different layers 
at different levels. The measurements were carried out in horizontal 
passages (about 2 m high) in the salt layers at two levels 660 and. 
540 m w.e. The work in the salt mine requires the apparatus to be 
built very thoroughly. Salt powder and moisture are deposited. during 
the work and therefore the whole apparatus must be tightly closed. 
These difficulties may be easier overcome if we use all-metal counters 
with earthed cathodes. Because of the very Jarge fluctuations of voltage 
in the mine, the apparatus had to be supplied with a stabilizer stabilizing: 
the voltage within wide limits. The magnetie stabilizer?) which has 
been used was very helpful in this matter. The special conditions of 
work also needed some additional arrangements to make work con- 
tinuous and safe (such as an automatie switeh tror too low voltage, 
a delayed switch for switching on the high tension for the register 
circuits, an electric clock set going by this delayed switch, eto.). 

The measurements of the ratio of threefold to double eoincidences 
were performed with an apparatus much improved in comparison 
with that of Barnóthy and Forró. The telescope consisted of three 
GM-counters 1, 2, 3 (Fig. 1). The middle counter 2 was of larger di- 


mensions than the two others. The telescope was protected against 
side showers by 6 anticoincidence counters. The double coincidences D 
of the counters 1, 3, and the threefold coincidences 7 of the counters 
1, 2, 3 were registered simultaneously by two separate P. O. registers. 
In this way we got a fairly small statistical error of the ratio T/D 
even for a not very great number of D coincidences. This is very impor- 
tant here because of the very low intensity of the cosmic radiation at 
the depths at which we have worked (the number of T coincidences 
at the lowest level is at the rate of about 1 per hour). In order to measure 
the double coincidences D the resolving time of our apparatus had 


z The authors are greatly indebted for the design and the construction of - 
this stabilizer to Mr. T. Czayka of the Cracow Mining Academy. 
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to be rather small. It was 2,6 usec, which we achieved in such a way 
that the rectangular pulse coming from the counter amplifier was 
differentiated with a small time constant in the circuit of the grid 
of the Rossi valves (Miesowiez 1947). The GM-counters we used were 
of the all-metal type ıilled with the usual argon-alcohol mixture. 
The dimensions of the counters were 


Counters 1 and 2. . . „43x65 cm}, 


Counter 3 and anticounters 5,2 x 70 T po tube a mn ihik 


The number of pulses per min. of these counters in the laboratory 
were about 650 for smaller and 900 for larger counters. The number 
of pulses per min. at the lowest level in the mine were about 70 and 
100 respectively. The low background rates of the counters were very 
convenient for our measurements. After protecting the counters 1 and 2 
by 5cm Pb at the lowest level, the background dropped to 16 per min. 


Measurements of the ratio of threefold to double and 
fourfold to threefold coincidences 


In our preliminary measurements performed at two depths 
(660 m w.e. and 540m w.e.) we found more double than threefold 
coincidences, proving in this way qualitatively even for those depths 
the existence of the effect found by Barnöthy and Forró at 1000 m w.e. 
We did not find however any difference between the number of three- 
fold and fourfold coincidences. The results of the measurements are 
given in Tables I, II, and III. In all tables the numbers of D-coinci- 
dences are given after subtracting the accidental coincidences. 


TABLE I 


Vertical telescope at a depth of 660 m e. w. 


number of coin- 
cidences and time 


average rate per - 
hour 


Threefold coince. T 321 1,34+ 0,05 
Double coince. D 443 1,86 + 0,06 s. 
Time (in hours) 239,46 


E 
57 


Ratio of threefold to 


double coine (72,5 2,7) Vp 
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TABLE II 


Vertical telescope at a depth of 540 m w.e. 


number of coin- 
cidences and time 


average Tate per 
hour 


Therefold coine. T 
Double coine. D 


2,08 + 0,07 
3,01 + 0,08 
Time (in hours) 


Ratio of threefold to 


- 0 
‘double coince. = (69,0+ 2,2) °/ 


SER 


TABLE III 


Vertical telescope at a depth of 540 m w. e. 


number of coinc. 
and time 


average rate per 
hour 


Fourfold coince. F 
Threefold coine. T 
Time (in hours) 


2,094 0,07 | 
2,104 0,08 : ; 3 


Ratio of fourfold to - 
threefold coinc. 
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ponent and A the efficiency of the GM-counter for this component. 
From our measurements of the ratio of F/T it follows that ANM, 
in other words the number of threefold coincidences produced by the 
non-ionizing component is negligible in comparison with the number 
of threefold coincidences produced by the ionizing component. Then, 
we may consider 7 as the intensity J of the ionizing component and 
write 

TIERE HET: 

Dc BN +2 ALT 


where A gives the number of double coincidences due to the non- 
ionizing component. Table IV gives a comparison of the values of 
A and T for two depths. 


TABLE IV 


coincidences due 
to ionizing com- 
ponent T 


eoincidences due 
to non-ionizing 
component A 


540 m w.e. 0,93/hour 
660 m w. e. 0,525, 


2,08/hour 
1,34 - ,, 


If A is considered as a measure of the intensity of the non- 
ionizing component of cosmic radiation, we should have a strong 
dependence on the depth and a high relative intensity. The ionizing 
component, however, shows that the intensity decreases with depth — 
in accordance with the measurements of other authors. For the expo- 
nent in the absorption formula we get from our figures the value 2,1, 
which taking into aecount our rather poor evaluation of the water 
equivalent, is quite consistent, within the limits of experimental errors, 
with the values given by other authors. ) 

The other characteristics of the component A correspond in 
general to those found by Barnóthy and Forró for radiation recorded 
by double coincidences. For the investigation of the angular depend- 
ence of the component A, measurements of threefold and double 
coincidences in a horizontal position of the telescope have been carried 
out. The results are given in Table V. 

From the figures given in it we can see that for a horizontal 
position of the telescope.the number of threefold coincidences is 
negligible and the number of double coincidences has, within the 
limits of experimental errors, the same value as A. This is an evidence 
for the isotropic character of the component A. 


| 
| 
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TABLE V 


The telescope in a horizontal position at 540 m w. e. 


number of coinci- 


i rate per hour 
dences and time 


Threefold coincid. T 
Double coincid. D 
Time (in hours) 


We have investigated also the influence of lead put between the 
counters on the ratio of T/D. The figures are given in Table VI. 


TABLE VI 


The telescope in a vertical position at 540 m w. e. 


'Thiekness of lead 
in mm 


From these figures we can see that the component A is absorbed 
with an absorption coefficient of about 1,4 cm! Pb. 


Efficiency of the counters for the component 4 


Barnóthy and Forró who evaluated the ratio T/D at 1000 m w. e. 
and found the value 0,05 interpret this figure as the efficiency of the 
counters for the radiation giving the coincidences. We think that 
this would be correet only if the telescope had not registered also 
a number of ionizing particles for which the efficieney of the GM-coun- 
ters is equal to 1,0, and for which the numbers of double and threefold 
coincidences are the same. For the efficiency A of the counter for 
radiation causing the excess of the double coincidences, we have A<T/D. 
On the other hand, in these circumstances we cannot at all evaluate 
the specific ionization as it was done by Barnóthy and Forró, who 
gave the value of 0,04 ions/cm N. T. P. air for this radiation. It is 
not known a priori whether this radiation reacts on the counter by 
producing ions in the gas of the counter, or whether it gives secondaries 
in the walls of the counter which initiate the discharge. In the first 
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case, the efficiency of the counter would depend on the specific 
ionization of the gas contained in it, but in the second case it would 
depend on the probability of emitting secondaries in the walls of the 
counter. 

For getting at least a rough evaluation of the efficiency of the 


' GM-counters for the radiation A, we compared the number of counts 


of a single counter caused by the radiation which is absorbed by 5 cm Pb 
with the number of double coincidences A. At 660 m w. e., for a counter 
of a length of 1=65 cm and a diameter of a=4,3 cm, the number 
of counts without screening was 80/min., with a screen of 5 cm Pb 
from all sides it was 16/min., thus the number of counts due to the 
radiation absorbed in 5 em Pb was 64/min. If we assume that this 
radiation is isotropic, we have a dependence between the number N 
of counts, the flux j of the radiation and the efficiency A of the counters 
for this radiation of the form ; 


N — nal (1 Hs a) Aj 


We get from it Aj — 0,022/sterad. cm? min. 
Let us now evaluate the probability of double coincidences due 
to this radiation, without any .discussion of the particulars of the 


mechanism of this process. Let A’ be the efficiency of the telescope 


. coincidences 


for double coincidences for this radiation, i. e., A’ is the probability 
that the particle or photon which discharges the first counter enters 
the second one and discharges it also. Assuming the radiation to be 
isotropie and paying due regard to the solid angle determined by the 
telescope, we get by simple integration for the number of double 


A= 24 jah are tang : 1 


where h—18 cm is the mutual distance of the counters in the double 


^. coincidence. Milencope: From this deb ee we get 


— 0,002, 


a value 20 times smaller than that given by Barnóthy and Forró. 


TNT, 


Conelusions. 


The results of this work therefore show the following properties 
of the radiation giving the excess of double over threefold coineidences: 

1. It is scarcely ^onizing. = 

2. It is isotropic. 

E» It is a by about 0, 5 cm Pb to 50 "Jo 
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All these properties seem to prove the faet that this radiation 
is of photon character, and the absorption coefficient as well as the 
efficieney of the eounters for this radiation also suggest that we may 
have here to do with some y-radiation of about 1 MeV energy, originating 
perhaps from naturai radioactive substances. 

It remains to explain the mechanism of the coincidences given 
by such a radiation. It is not improbable that we are confronted here 
with an example of some compound Compton-effect. As a matter 


of fact, it might be possible that a photon entering the first counter 


and discharging it might enter afterwards the second counter and 
give off there another electron which discharges this second counter. 
It may be seen from our evaluation of the efficiency of the counters 


for that radiation — assuming that the discharge in the counter is : 
started by photons — that we get thus quite a reasonable figure for © 


the efficiency for the double coincidences. Of course, in this case 4'« 4, 
because the loss of energy of the photons in the first counter causes 
a drop in the efficiency and the photons are dispersed producing in 
this way also an additional inefficiency. 


We wish to thank the Science Department of the Ministry of | 
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PROPAGATION OF A CUT-OFF TRAIN OF DE 
BROGLIE WAVES 


By A. RUBINOWICZ, Institute of Theoretical Mechanics, University 
of Warsaw, Warsaw 


(received July 17, 1949) 


The propagation of a cut-off train of monochromatie de Broglie waves has 
been investigated. Formulae have been developed describing completely the arrival 
of the wave front, the forerunner, and the main wave. 


In the present paper we deal with the following problem of the 
propagation of a cut-off train of plane de Broglie waves. At the time 
1—0 we have at the negative x-axis a train of monochromatic de Broglie 
waves of frequency »,, that is 

Zh (e-&) (1) 
with a phase velocity u, given by 


C 


where »,=m,¢?/h is the frequency corresponding to the de Broglie 
wave length of a particle of rest mass m,; a is an arbitrary phase 
constant. Our problem is to determine the propagation of these waves 
for t>0 and >Q. 

This problem is formally a special case of Sommerfeld’s (1914, 
cf. also L. Brillouin 1914) problem of the propagation of an electro- 
magnetic cut-off wave train in a dispersive medium, namely the one 
dealing with free electrons without damping. This specialization enables 
us to obtain convergent series developments of the exact solution, 
whereas in Sommerfeld's general case it was only possible to discuss 
it approximately. 


(2) 


§ 1. Solution of the problem and its series develop- 
# 
ments. The integral l 
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taken in the upper half of the complex »-plane along the path (U), 
fig. 1, from +00 to —oo of the real axis defines the function 


0 for £« 0, 
«tL. g- nit for 10. 


For t<0 we can namely deform the path (U) into the infinity of the 
upper half of the »-plane where exp (—2zi»t) for t< 0 vanishes exponen- 
tially. For t>0, however, we can deform (U) into the infinity of the 
lower half of the complex »-plane, but in doing so, we must pass the 
pole at v=»,, which gives a residue. 


v- plane 


Fig. 1. 


To a corpuscular wave of frequency » corresponds a phase 
velocity u given by 


e 1-3)" (3) 


We expect, therefore, the solution of our problem to be 
ela dy 
Í (t, x) = Oni e 
(U) 


(4) 


y A Vo : 
where 


w= — 2niv (:—2) 22 (n—5 er). (s 


Owing to the fact that w in (3) is positive, we must give to the 
square root in (5) a positive sign for real »-values v >». The points 
?»—-k», in (5) are branch points. In order to have to do with a definite 
branch of w, we place a branch cut on the real äxis of the »-plane 
between these two points. 


First of all, we can show that f((,z) vanishes for i<= , that is 


for B>1, if 


p--—. (6) 
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This means that the velocity of the wave front cannot be greater than 
the velocity of light. In the case f >1 we can indeed deform the path 
of integration (U) into the infinity of the upper half-plane, because 


for large »-values the function w, eq. (5), is given by w= — 2riv (- 2) 


Hence, 
f(t, æ) Oe Ion > 


For 8—1, however, we can deform (U) into the infinity of the 
lower half-plane, where, in this case, w vanisbes. But then there remains 
still a path of integration (U,) surrounding the branch cut between 
—7?, and +», and the pole at v—»,. For B —1 we can therefore use instead 
of (U) the path of integration (U,). 

The meaning of f is as follows: according to (6), the value 
of cß=x/t determines the velocity of a particle which would reach 
the point x at the time t starting from x—0 at the time £— 0. 

In order to obtain convergent developments for f(t,z) in the case 
of B<1, we deform (U,) into a path of integration (E), fig. 2, on which 


Fig. 2. 


the real part of w, eq. (5), vanishes. (E) is an ellipse with its centre 
at v=0 and with the principal semi-axes a 


"n depu s m di (7) 
per ar NP 


a and b are situated on the real and imaginary axes of the complex 
v-plane respectively. (Æ) can therefore be represented by 


a -— 


Vi oe 


dice [à 


yz (cos y + i sin y) (0x v « 2a). 
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Because on (E) 


x 1—P? 
w= w, =— iz cos y, where ne! ^. (8) 
we have 
pui 
is — sin ip eos 
f(t,o) — = an dy. (9) 
Tg cos p+ é8 sin p——2 iR 
1 
Putting 
V Yo (22) 
— nett 2£2]/ IE EE Mo ern 10 
yet y Tap Yo eV v , (10) 
we obtain 
— sin y +76 cos y = @—1 
cos y+ if sin y — + V1— p aztia 
1 : 1 
Yo 1 1 (11) 
=i/1+-%-+— 
6 — yo Yo si) 
Yo, 
so that finally 
ar 
ei® Yo 1 1 , : 
1,2) —.— pe"|l---——--— d residue). 12 
fta) ad tpe teme. — 02 
yo 


„+ (residue)“ means here that the residue of the integral (4) which 
is given by the eqs. (1) and (2) must be added to the integral in (12) 
in the case of the pole at v—», being situated outside the path (E). 

As the frequencies v, and », are determined by hv— mc and 
hv,— msc? respectively, where m, is the rest mass and m the mass of 
the moving particle, v, is always greater than »,, so that the pole 
at v= lies outside of the branch cut. For ß—1, according to (7), the 
semiaxes of (E) become infinite, and for ß—0 the path (E) approaches 
the branch cut. If », is equal to the greater semiaxis a the path of 
integration (E) passes through the pole at v—»,. Denoting the corre- 
sponding value of f by f,, we get, according to (7), 


^ 
V = SS. 13 
Inserting this into the eq. (10) for yọ, we get 


n= |/ith. 114) 
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Using (12), we can now obtain representations for f(t,~) by means 
of convergent series of Bessel functions 


27 


jn 
ACE 4 a; [ mtem icis f nmmennis (15) 


0 


We deal with the case 1 >£ >, first. According to (10) and (14), 
we have then +co>y>yp, so that owing to y—Jc| 


Too >|] >ro- 


The path of intergration (E) contains in this case the pole at v—», in 
its inside, so that there does not appear any residue in (12). According 
to »,>», and (10) we have y,>1. Under the supposition that |] >yo, 


1 
we get have also |¢ Pay , 80 that the bracket expression under the 
0 


sign of integration in (12) can be expanded into the power-series 


o1 
PETIT ty. 


n—i 


Hence, with the help of (12) and (15) we get for the function f(t,«) 
the expression 


f(t, æ) = et^ (J a(z) + 2 (—i) J ale) y" (yg- Yo 7))- (16) 
The moment t, given by ß,=x/et, defines the velocity v9 — x/tp= cf, . 


2] 
As, according to (13) there is f= (1 -(&j) , we have 


07 
y, \2\"/a 
«e - CI) 


=. 


and hence, using (2), 


ty is therefore the group velocity corresponding to the phase velocity ug. 
At the time f£, there arrives at the point c that part of the wave train 
which travels with the group) velocity Vo, i. e. with the velocity of the 
incident particles. 

This means that the development (16) describes the events for 
a given point # from;the arrival of the wave front to the arrival of 
the waves travelling with the particle velocity v,. This development 
can be used also for times t greater than tọ, that is for fa >f 2 0. But 
in this case the convergence is very slow. 
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For these f-values we get, however, a good convergent series 
Yo 1 

and — ———— 
— Yo vo 6— (1/70) 
power-series in C/y, and 1/&y, respectively. Taking into account 
that now the pole at v=» lies outside of the path of integration (E) 
and therefore the residue must be added, we obtain 


by developing in (12) the expressions f into 


TUIS TAN [mese - Sgr}. an 


n—i 


(17) also can be used for all f-values in the interval (0,1). From 
the equality of (16) and (17) we obtain a development of the incident 
wave given by (1) and (2) in a series of Bessel functions 


mme se G1) + Ste ((22)'+(2)') 


Yo 


This expansion follows also directly from the generating function 
of Bessel functions (Watson 1944) 


co 


Jo(2)+ Inte) ((4- (—t)) 


n=1 


ee _ 


if we put here t= —iyo/y. 

For large values of 2 we need an asymptotic expression for the 
wave function f(t,z). We can get it by applying to (9) the principle of 
stationary phase. For large z-values the oscillations of e% are very 
rapid in comparison with the changes of the remaining part of the 
integrand in (9). The phase wọ has stationary values for y=0 and y= x. 


2 2 
Using the approximations cosy=1— + and cos (z+ 9») =—1+4 + 


in the neighbourhood of these points, we get in the usual way 


: ia 12-5) +i 2-7) 
f(t, 9) = enr cup aut cialis ti + (residue). (18) 


Vax Ve \»—N-# | 4 +%VI—P 


Formula (18) can be used only if the path of integration (E) does 
not pass through the neighbourhood of the pole at v, independently of 
the fact whether », lies inside or outside of (E). In the second case, 
that is for f,>ß>0, we must, however, take into account the 
residue at v=»,. 

Asymptotic expansions containing descendig powers of e and the 
expression (18) as the first term can be obtained by developing the 
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; 2 2 
product of exp | —i2 (eos v(t -£J) Or exp E i2 (cos (x4-9)4- i-£J 
and the expression (11) in a series of ascending powers of y and g re- 
spectively (an der Corput 1939). We can also get the expression (18) by 


the method of steepest descent. 


$ 2. Discussion of the solution. In order to have an idea 
about the order of magnitude of z we introduce in (8) by 


h 4 
dere (19) 
TC 9 


the Compton wave length A of the incident particles. Thus we get 


æ y1—pg 

2=27 ina? IER (20) 
For electrons, for example, is A=2,4 10-10 cm, and therefore is z 
a large number for z4-0 and f<1. 

First of all, we consider the arrival of the wave front at a given 
point x. It comes at the instant determined by ß=1, and hence 
by 2—0 and y=oo. As J,(z) is the only Bessel function which does not 
disappear for 2— 0, it follows from (16) that for f—1 


f(t, ©) = eJ (2). (21) 


The intensity |f(t,z)|? of the arriving wave front is therefore of the 
order of magnitude of 1, that is of that of the incident wave. More- 
over, (21) is independent of the frequency of the incident wave. We are 
-hence compelled to assume that the events at the wave front are deter- 
mined only by the process of the cutting off of the incident wave. 

We ask now whether it is possible to assign a wave lenght 4 and 
a frequency n to the front wave (21). The distances of two consecutive 
zero points of the Bessel functions J„(2) are of the order of magnitude 
of x. The variable z can therefore play the role of a phase. If e would 


be equal to 2z (7 -i) then the wave length A and the frequency 
n=1/T would be given by 2: d 
| 19 08 1) d: ee | 
Dana ee Ts Dean: e 


But following L. Brillouin (1914) we can use (22) as a formal defi- 
nition of A and n respectively. Using f-—«w[ct we get then from 
(20) and (22) 

(ur d Ay ne (23) 
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Ais the de Broglie wave length of particles moving with the velocity cf, 
n is the frequency associated with A by the de Broglie phase-velocity c/p. 
Owing to (19), we have An=An»/ß=e/ß indeed. At the moment of 
the arrival of the wave front (f —1) are 4— 0 and n=oo; afterwards A 
increases and n decreases. 

The moment the wave front has passed the point x (that is for 
B<1) the value of z increases so much that we can use the asymptotic 
expression (18), first for 8,<ß<1 without the residue. This part 
of the wave f(t,z) is called the forerunner or precursor. Owing to Vz 
in the denominator the function f(t,x) decreases until f comes into 
the neighbourhood of f. Then, according to (13), the resonance denomi- 
nator of the first term in brackets becomes very small, so that f(t,) 
increases again. But at the same time the asymptotic expression (18) 
ceases to be applicable, because the path of integration (E) passes 
through the neighbourhood of the pole at »—»,. In this moment the 
„main wave“ given by the residue arrives. Subsequently, the contribution 
of the first term in the brackets, which at v—», has changed its sign, 
becomes less important. With the advance of time (f—0) the value | 
of Vz approaches infinity, so that practically only the main wave 
remains. 

The expressions (23) for A and « are valid also for the first term 
in the asymptotic approximation (18). With decreasing f, the. wave 
length A increases, until at the moment (f= f$) of the arrival of the 
main wave it becomes equal to the de Broglie wave length of the 
main wave. After this möment A inereases further. 
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COMPARISON METHOD OF MEASURING THE ANGU- 
LAR CORRELATION OF GAMMA-RAYS 


By Jerzy GIERULA, Physical Laboratory of the Mining Academy, 
Kraków 


(received October 30, 1949) 


From theoretical considerations (Hamilton 1940, Goertzel 1946) 
it follows that when two y-quanta are emitted in cascade in the 
proeess of nuclear de-excitation the relative probability that the 
second quantum will be emitted at an angle O with the first is 


l 
W(0)—1 +> A, cos?! 0, 
=i 3 


where 22 is the highest multipole order occuring. 

Up to now, the investigations of the angular correlation of y-rays 
(Beringer 1943, Boehm 1949, Brady 1947, Good 1946, Kikuchi 
1942) were carried out by observing the coincidences between a coun- 
ter A, fixed with respect to the sample P, and a movable counter B 
(Fig. 1). As these measurements extend as a rule over 
long periods of time, they require a far reaching stability ©) 
of the working conditions of the counters and an exact , 
reproducibility of the geometrical configuration. If the ir 
effects of the angular correlätion are small (as is to ^ H 
be expected in the majority of cases) a very slight D 
instability of the experimental conditions may spoil the ar 
measurements to such an extent as to make them wholly Fig. 1 
unreliable. 

In this note a comparison method of determining the angular 
correlation is suggested which is more or less insensitive to the insta- 
bility of the working conditions of the counters and of their spatial 
disposition. During the whole course of a series of measurements two 
countersA and B (Fig/2) are fixed, while a third one C is placed in 
position 1 during the first part of the series and afterwards in position 2. 
The first part of the measurements consists in the simultaneous counting - 
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of the (A,C,)- and (B,C,)-coincidences (counter C in position 1). The 
numbers of these coincidences are given by 


(A, €) —N,0,0, £46; V ag W (909), (B, C) —N 05%. ££; V gc W (180°), 


where N, is the number of decays during the measurements, c are 
the respective solid angles, e the efficiencies of the respective counters, 
V are the respective losses of coincidences, and W(@) is the relative 
probability of emission of two successive quanta at an angle ©. The 
ratio of the above numbers is given by 
(A, Cy o 8,77, W902) a) 
(B, Cara W (180°) 
Similarly, the second part of the measurements with counter C 


in position 2 gives 


(A, 0,)=N,0,@ 0, E48cV 4c W (1809), (B, C,) — N,0 p00, £g&c V ge W (90°), 
and their ratio 
(AS 03) aren ace WiCESOR) , 
(B,O) "eye, V, WOO) ` (2) 


From (1) and (2) we get 
VE /4,0) _ W(90°) (3) 


(B,O)/ (B,C) W(809' 

a quantity which is independent in a high degree of the fortuitously 

changing working conditions of the whole coincidence apparatus. 
A In Fig. 2, the counters are shown in position 
C for a measurement of the ratio of the probabilities 
E : < at 90° and 180°, but it is easily seen that this ratio 
oo. may be measured for any other two angles by an 
. | appropriate change in the angular setting of the 


ons. 


oy counters A, B, and C. 
"ice Of course, the same method may be also used 
Fig. 2 for angular correlation measurements of (ß,y)- and 


(8, B)-processes. The errors introduced by the absorption of the f-particles 
in the sample may be avoided by a steady rotation of the sample 
around its axis. 
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DIPOLE MOMENT OF CARBONYL CYANIDE 


By M. PUCHALIK, Institute of Physics of the L. Waryński Silesian 
Academy of Medicine, Rokitnica 


(Received March 16, 1950) 


The dipole moment of carbonyl cyanide has been investigated. This com- 
pound, for the first time obtained by Prof. Malachowski and collaborators (1937) 
has an analogous structure to phosgene. It has been ascertained that in carbonyl 
cyanide, just as in phosgene, there appears a compensation of dipole moments of 
two strongly polar groups. 


The investigations of the polar properties of phosgene by Smyth 
and Me Alpine (1934) and Le Févre-and Le Févre (1935) led to 
the following conclusions: (1) The molecule of phosgene forms a dipole 
with an electric moment— 1,18 D (1934) or 1,099 D (1935). (2) The 
dipole moment of COC], may be calculated from the values of the 
moments of the two groups — CO and -C-Cl. (3) The relatively low 
dipole moment of phosgene in spite of the presence of strongly polar 
groups is due to mutual compensation of both their moments. 

It seemed interesting to investigate the polar properties of car- 
bonyl cyanide, as this compound has an analogous chemical structure 
to phosgene. Carbonyl cyanide was synthesized for the first time by 
Prof. R. Malachowski and coll. (1937). This remarkable substance 
shows great chemical activity, e. g., it reacts violently with water 
giving CO, and HCN. As hydrogen cyanide shows a rather high dipole 
moment (u= 2,5 D) all precautions were to be taken to avoid humidity 
of the substance as well as of the apparatus. Even small amounts of 
water would strongly influence the results of the dipole-moment 
measurements. A special capacitor was used which made it possible 
to protect the examined substance from any influence of humidity. 


The construction of this capacitor is shown schematically in fig. 1. 


pae 


In principle, it is a cäpacitor with two dielectrics: glass and the in- 
vestigated substance. Before use the capacitor was always carefully 
dried with a current of air passing through H,SO,, CaCl, and P,0,. 
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The capacity of such a capacitor when connected to the measuring 
circuit is given by the formula 
Qu dM (1) 
Be+D’ 
where A, Band D are constants (A=BD), C, is the capacity of the 
connecting wires, and e the dielectrie constant of the examined sub- 


stance. C, is very small in comparison with C. From (1) we get 
__D(C=6,) 
STMZERIO 0): 
Let e and ej, denote the dielectric constants of the solvent and 
the dilute solution respectively. As the difference 4C of the capacities 
for e=&,, and e=s, is small, we may write 


(2) 


Be 
Er —-& = (0,6 AC 
or, neglecting higher powers of C)/C, 
Be 
£5 — E+ AC 


The dielectric constants of the two non-polar solvents used 
(CCl, and C,H, &;5:— 2,245 and 2,278 respectively) differed only slightly 


Fig. 1. 


and, therefore, the bracketed expressions in the denominator may 
be considered as a constant, say H/B, yielding finally 


e | 
f= +H; 40, (3) 
1 


where H is a constant to be determined by experiment. This was done by 
calibrating the capacitor with benzene, carbon tetrachloride, and hexane. 


u 
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The determination of the dipole moment of CO(CN), was performed 
by the method of dilute solutions; those of the dielectric constant 
by means of the bridge method at audio frequency. A general scheme 
of the experimental arrangement may be seen from fig. 2. 


4V 
(00V 


Fig. 2. 


The results of the measurements are shown in the following 
tables t), where f, is the mulecular concentration of the carbonyl 
cyanide, e, and d,, are respectively the dielectric constant and the density 
of the solution, and P, is the molecular polarization of carbonyl cyanide. 


Table I 
CO(CN), in C,H, T=25°C, R=16,93 em? 


dy, g/cm? 


0,8927 


0,8823 


0,8768 


0,8756 


F 
PP?-—64cm?, w=1,5D45%, 


1 Preliminary results were communicated at the IX Meeting of Polish Phys- 
icists at Wilno, October 1938. 
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Table II 
CO(CN) in CCl, 2=25°C 


P>=64 cm, a o DESSY 


Thus, the analogy of the chemical structure of COC], and CO(CN), 
reveals itself also in the similarity of their polar properties. Like 
phosgene, carbonyl cyanide shows in spite of the presence of two 
strongly polar groups a comparatively low value of its dipole moment. 
An approximate calculation for carbonyl cyanide based on the known 
group-moments for— CO and =C-CN gives a value between 1 D 
and 2D. 

The only difference in the polar character of these substances 
consists in the fact that the moment of the carbonyl group is greater 
in the case of phosgene and smaller in the case of carbonyl cyanide 
than the geometrical sum of the moments of the C-Cl snd the C-CN 
groups respectively. - 

These investigations were carried out in the Institute of Experi- 
mental Physies of the University, Lwów, in 1938—1939. By comparison 
of their results with modern litterature the author has got the impression 
that their value remains unchanged. 

The author wishes to express his gratitude to Professor S. Loria 
for his most valuable advice in the course of this work. 
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DIPOLE MOMENTS AND FREE ROTATION IN 
MOLECULES 


By M. PUCHALIK, Institute of Physics of the L. Waryński Silesian 
Academy of Medicine, Rokitnica 


(Received March 16, 1950). 


To test some consequences of the concept of free inner rotation in molecules 
(a) the mean square dipole-moments of some members of the homologous series 
C;H,(CH.)nCl, C,H,(CH,)„Cl and C,H4,,(CH;)4Cl, and (b) the temperature dependence 
of the mean square dipole-moment of benzyl chloride have been investigated. The 
results seem to confirm the theoretical predictions of Zahn (1932) and the results 
of some, not numerous, experiments of other authors. 


The concept of free rotation in molecules is based on the fact 
that in many cases where -the structure formulas provide for the 
existence of geometrical isomerism, such an isomerism does not appear. 
It is obvious that in polar molecules the freedom of rotation cannot 
be complete; the mutual repulsion of the polar groups makes the trans 
configuration of the molecule the most probable one, as this config- 
uration is connected with a minimum of potential energy. According 
to known chemical facts one must assume that the difference between 
the minimum and the maximum potential energy is small and the 
energy of the molecular motion is sufficient to be able to change the 
trans configuration into the cis configuration and vice versa. i 

It may be inferred on theoretical grounds that the mean square 
dipole-moment of molecules which form a stable system with an 
axis of completely free rotation is given by the following simple 
formula (Zahn 1932) 


TA E 


where u, is the geometrical sum of the rigid group- -moments and the 
components of the rotating group-moments in the direction of the 
axis of free rotation. It is evident that during the rotation these compo- 
nents keep their magnitude and their direction constant. u, are the 
components of the rotating polar groups perpendieular to the axis 
of rotation. 
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In the general case in which the single axis of rotation can rotate 
freely, the value of u2 may be calculated by successive application of 
formula (1). 

The determination of the degree of the freedom of rotation may 
be of great assistance in solving different problems connected with 
the chemical structure of molecules. It may be performed by means 
of the following considerations: The values of the mean square dipole- 
moment of a chemical compound with an incompletely free inner 
rotation depends on the temperature as well as on the mutual distance 
of the rotating polar groups. The increase of both these factors makes 
the freedom of rotation more and more complete and at the same 
time the mean square dipole-moment tends towards a limiting value 
given by formula (1). For rigid molecules and for molecules with 
a complete freedom of rotation the mean square dipole-moment is 
independent of the temperature. 

In the ease of a homologous series having polar groups at both 
ends of a carbon chain the lengthening of the carbon chain increases 
the distance of the polar groups, so that the higher members of such 
a series ought to show nearly constant values of their dipole moments, 
which corresponds to a practically completely free rotation. This was 
confirmed by the investigation of Riedinger (1938) on the members 
of the homologous series Br(CH;),Br and C,H;(CH;),C,H;. 

More extensive systematic researches on this matter seemed to 
be advisable, and so the author decided to verify the validity of the 
above mentioned considerations on the inner rotations in molecules 
by experimental investigation of some members of other homologous 
series. The determination of už were carried out by the method of 
dilute solutions, benzene and carbon tetrachloride being used as non- 
polar solvents. The dielectrie constant was measured by. means of 
a bridge method described in a previous paper (Puchalik 1950). 

The members of the following homologous series were investi- 
gated: (1) C;H,(CH,),Cl, (n=0,1,2), (2) C,H,(CH,)„Cl (n=0,1,2), 
(3) CH (CH;)401 (n— 0). i 

The mean square dipole-moments of the following chemical 
compounds belonging to the above series were determined for the 

first time: (1) Cyelo-penthyl chloride, cyclopenthyl-methyl chloride, 
eyclopenthyl-ethyl chloride, (2) -phenyl-ethyl chloride; for chloro- 
benzene and benzyl chloride the results obtained are in good agreement 
with previous determinations of other authors, (3) Cyclohexyl chloride; 


the mean square dipole-moment of 2-p-cymenomethyl chloride has 
also been measured. 


Dipole Moments and Free Rotation in Molecules 95 


The results of the measurements are shown in the following four 
tables, where f, is the molecular concentration of the dissolved polar 
substance and R is the molecular refraction. 


Table I 
Homologous series C,H,(CH,)„C1 
A. Cyelopenthyl chloride in CCl, 
t= 25° CP Ny=1,4488 d=0,99458 g/em3 R= 28,15 cm? 
di, 


fz E12 


1,5635 


115 


V= 1,98 D + 59/, 


B. Cyclopenthyl-methyl chloride in CCl, 
t= 25°C, R= 32,78-0m3 


ds g/cm? 


2,749 1,5297 


2,438 1,5585 
2,377 1,5718 


V2 = 2,35 D+ 5%), 


C. Cyclopenthyl-ethyl chloride in CC], : 
t=25°C, R= 37,39 cm? 


V5 = 2,38 D+ 5%, 


96 M. Puchalik 


Table II 
Homologous series C,H;(CH;)4Cl 
A. Benzene chloride in benzene 


t=25°C, np= 1,5219, d= 1,0952 g/cm?, R= 31,26 em? 


E12 
3,163 
2,838 0,9138 


2,649 0,8993 
2,363 0,8875 


P2=84cm?, V= 1,58 D 4: 5/, 


0,1593 


0,0647 


B. Benzyl chloride in benzene 


t=25°C, n,—1,5362, d= 1,09259 g/cm?, ` R= 36,12 cm? 


12 
3,374 0,9408 


Vis — 1,86 D + 59, 


C. #-phenyl-ethyl chloride in CCl, 
(—25 C &=5,43 nj—1,9280 d— 1,0655 g/em* R= 37,73 cm? 


Surface tension ø= 42,7 dyn/em Viscosity u= 0,0192 cm! g/seki 


0,2310 


0,1345 


0,0742 


0,0449 


V2 = 1,92 D + 50], 
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Table III 
Cyclohexyl chloride CCl, 
t=25°C, n,=1,4606, d=0,9910 g/cm? R= 32,80 cm? 


Vut = 2,2 D45 


0,0973 


12 
0,0509 5 


Table IV 
Karvakryl chloride in CCl, 
t=25°C, np= 1,5233, d— 1,0104 g/cm? R= 51,29 cm? 


dia g/cm? 
> 
, 


Dum pre 


V®=2,1 D+50), 


12 
, 


These experimental results show that: (a) In good agreement 
with the above mentioned theoretical considerations, the change of 
the mean square dipole-moment for the two first homologous series 
becomes small with increasing n, which means that the freedom of 
rotation becomes soon practically complete. For the second series 
the limiting value of the mean square dipole-moment can be calculated 
and one gets the value 1,9 D, possessed already by f-phenyLethyl- — 
chloride. (b) The mean square dipole-moment of benzyl chloride is 
less than of cyclohexyl chloride and cyclopenthyl chloride. The change 
of the benzene ring to an alicyclic ring produces thus an increase of 
the mean square dipole-moment. Unfortunately, the number of the 
investigated members of each series was rather small. 

These investigations were performed in the Institute of Experi- 
mental Physies of the University, Lwów, but were interrupted in 1941 
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by the outbreak of war between the Soviet Union and Germany. By 
comparing their results with modern literature, it has been ascertained 
that they have retained their value up to now, as the state of the 
experimental knowledge of the dependence of the mean square dipole- 
moments on the freedom of rotation seems to have advanced but 
little during these last years. 

It would be interesting to be able to verify the above mentioned 
conclusions by means of measurements of further members of the 
homologous series investigated in Lwöw; unfortunately, it is now 
very difficult to obtain the catalyzers needed for the synthesis of 
these compounds. 

As the increase of the temperature has a similar effect upon 
the mean square dipole-moment of the molecules with an incomplete 
freedom of rotation as the increase of the carbon-chain length, the 
author has decided to investigate the dependence on temperature 
of the mean square dipole-moment of benzyl chloride. These researches 
have been carried out in the Institute of Medical Physics of the Silesian 
Medieal Academy. The determination of the mean square dipole- 
moment in the temperature intervall from 20? C to 80? C have been 
performed by the method of dilute solutions, refined petroleum being 
used as non-polar solvent. The polarization of benzyl chloride has 
been calculated by means of the formula (Müller 1933) 


1 &o—1 &— 1 
2 


Fin, Sad "esi (2) 
where n, is the number of gram molecules of the dissolved polar 
substance, and V,, are respectively the dielectrie constant and the 
volume of the solution, and V, are the same quantities for the non- 
polar solvent. _ ; 

The dielectric constant has been measured by means of the 
bridge method, a measuring capacitor being used as in the previous 
work of the author (Puchalik 1950). The measuring apparatus pro- 
duced by the establishments „Duo“ at Gliwice is similar to the one 
previously used. Its schematic diagram is shown in Fig. 1. Modern 
vacuum tubes and feeding from alternating current network has been 
used. The accuracy of the determination of the mean square dipole- 
moment was about 3 per cent. The whole arrangement was tested 
by means of a solution of benzyl chloride in benzene. 


The results of the measurements are shown in the following 
tables. 
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Eh, 


N26 0- 524 


ee = (3 
E sr 


Fig. 1. 
Table V. t=20°C 


V2 = 1,69 D 
Table VI. t=40°C 


99 


ne 


———" 
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Table VII. t=60°C 


EN 
| 
Ce oe p 


i z1,79 D 


0,1920 


Table VIII. t=80°C 


ABER TI 
82,86 2,360 13,02 
84,10 2,526 16,15 


0,1114 66,34 2,583 11,51 
76,01 2,114 13,55 
0,0290 1,994 15,73 


Vu? = 1,85 D 


These experimental results show that: 


(1) In-the case of petroleum as non-polar solvent one obtains 
for the mean square dipole-moment of benzyl chloride a somewhat 
smaller value than by using benzene as solvent. 


(2) The mean square dipole-moment of benzyl d varies 
very little with temperature. 


(3) Thus, one is drawn to the conclusion that the freedom of rota- 
tion in the benzyl chloride molecule is not yet complete. 

The author wishes to express his gratitude to Professor S. Loria, 
Head of the Department of Experimental Physics of the University, 
Lwöw, for his valuable advice during. the first part of this work, and 
to the late Professor S. Pilat, Head of the Institute of Potrdfenma 
Tolosam at Lwów for his valuable ci and for, supplying 
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him with all the needed substances. He takes also pleasure in thanking 
Professor W. Lesnianski, Head of the Department of Organie Technol- 
 ogy of the Silesian Institute of Technology, Gliwice, and Dr. W. Ma- 
_ronski for supplying him with benzyl chloride. 
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LETTERS TO THE EDITOR 


The Board of Editors does not hold itself responsible for the opinions expressed by 
the correspondents 


On the Theory of Non-Local Fields 
JERZY RAYSKI 
Physical Institute, Nicholas Copernicus University. Torun 


April 30, 1950 


Recently Yukawa! initiated a relativistic theory of particles with finite radius 
based on Born’s ideas of non-localizability and reciprocity. The non-local field 
operator U satisfies the following commutation relations 


(du, [du ,U]] - x* U, [£us [xu ,U]] — 4? U, [du »[%p, U]]-— 0, ; (1) 
where d, is the displacement operator ] 
[zu ; dy] — 9n» - (2) 


In the representation where x, are diagonal U is a matrix Qz'|U|x'» which 
.may be regarded also as a function U(X,r) with 


Xp=1fAHe+ tu) Ya (Ya — Rn). (3) 
U(X,r) satisfies the equations 
di 2 , 
& i ar U(X,r)=9, (4^) 
na =, (r 32) = (4^) 
? , IX 


Equations (4") may be treated as supplementary conditions. The general 
solution of (4) is 


U(X,7) = T2 (u(r) eX + wk tet 70) 8 ((ker)) 5 ((r-7)— 2?) (5) 


with (k-k)=—x?. 

If the scalar U(X,r) shall describe spinless particles, then no direction can be 
priviledged except ky, so that uz(r), wg(r) may depend only on (r-r) and (r-k), which 
means that they do not depend on r, at all, due to the supplementary conditions. 
Thus, for spin zero we may write in (5) uz, wy instead of ur(r), wr(r). 

We may establish a correspondence between a local quantity a and a non- 


local A in the following way 
a(X) —» f dir qa! | Ala”. (6) 
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If A itself is local, then it is of the form 
«a | A| > = ale’) 6(x’— 2x’) =a(X) é (7) 


and the correspondence becomes an identity. From (6), there is a correspondence 
between the traditional local field function (X) and the non-local U(X.r): 


(x) —> F dirU(X,r) (7) 
or 


» (ag ei: X) + bž e iX) =, (uj e X) + wie X) Sp, (7^) 
k k 


where S4 is the surface of the ellipsoid 


Sa f dr ölk) (ra. (8) 


(7’) enables us to write immediately the commutation relations for U(X,r). We 
see that u,;S; must satisfy the same commutation relations as ag. Hence, we may put 


UX- >) g, c4 bf MEX) (7) &((-7) — 22), (9) 
k 


where a;,b; are the well-known matrices satisfying the commutation relations 


*) *1 b 9, 
[ak a1 ]  [ba. br ] Ve V (10) 
while all other commutators vanish. 

In order to have a complete theory we must find the important bilinear expres- 
sions for the Hamiltonian energy density, charge density etc. To do so we have to 
replace the products of y, w* and their derivatives by corresponding expressions 
in terms of U and U*. The correspondence rule (6) yields 


y*(X) p(X) > i dirqz'|U*U|z/^» = f dir [A diz'"x'(U*|z"»«x'"|U|x^». (11) 


With the aid of the notations 


= å y? = poutas yr migi 
2 (12) 
r= x’ — r”, r = z'—qc'", t” = qus ee 
(11) becomes 
v*(X) v (X) > f dir f dax" U*(X',r')U(X" v") 
(13) 


i f dar’ if dr’ U*(X4 12,0’) U(X 4- v], v"). 


By rewritting the traditional Hamiltonian for a scalar local field by means 
of (13) and using (9) and (10), we easily find that the new Hamiltonian has the same 
eigenvalues as the former and yields the correct field equations (4) in the operator 
form. The traditional expression for the charge density, rewritten with the aid 
of (13), yields also, if integrated, the correct eigenvalues for the total charge. Thus, 


the theory of spinless non-local fields in vacuo yields the same results as the local - 
field formalism. 
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Interaction between two spinless fields. Let us consider the simplest 
possible type of interaction between a complex and a real scelar field: 


A(x) = gap*(x) p(x) y(x). (14) 


This interaction energy density is subject to the well-known inconsistencies 
inherent in any local field theory. To take into account the non-localizability we have 
again to apply the rule (6). We limit ourselves to two special cases: 

(1) g local, v non-local. In this case the correspondence rule (6) yields 


H(X)=9 far’ b a" USC er" (X E ) cer.n. (15) 


(2) v local, p non-local. In this case we find 
H4X)—g j dirq*(X—T7[2) V(X,7) (X 4-72), (16) 


where V(X,r) is the non-local neutral field. The interaction appears suitably spread 
out over the ellipsoids representing thé non-local particles. 
Les us ealeulate from (16) the self-energy of a neutral non-local particle at rest: 


Heey= > PP, (17) 


Here, the initial state, denoted by 0, represents one neutral particle at rest 
> 
and the intermediate state, denoted by p, represents a pair of charged particles 


> > 
with momenta p and —p. The matrix element Hs differs from the usual one by 
. & factor a: 


a= F dire PD NUR § (479) d ((r-r)— A?) (18) 
0 - 
with 
> > 
P,=(?|VeP +I), q,—C—» | VP+ 27. 


S, is now simply the surface of a sphere with radius 4. Performing the integration, 
we get 


. sin Ap 
Css rae 


Thus, instead of a divergent integral over dp we get a convergent one with the 
„damping factor“ a? in the integrand. 

We have shown that the non-local field may be treated by the usual Hamil- 
tonian method. Contrary to the opinion expressed by Yukawa! it is neither 
necessary to reject the Schrödinger equation not to replace the Hamiltonian forma- 
lism by the S-matrix scheme. We hope that the well-known convergence difficulties 
: may well be overcome in terms of this new formalism. Moreover, this theory may 
be able to put the problerh of elementary particles with higher spin in an entirely 
new light. 


(19) 


1 Yukawa H., Phys. Rev., 77, 219 (1950). 
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